mas deve fazer notar-se que o caso mais corrente é 
o da precipitação dum hidróxido. 
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Fig. 3 Representação de equilibrios 
quimicos dependentes do ph 


Na Fig. 3 representa-se a família de linhas corres- 
pondente ao equilíbrio Fe(OH),/Fe”*. Estas linhas 
dividem o plano E-pH em duas regiões: a de estabi- 
lidade termodinâmica dos iões Fe” e a de estabilidade 
termodinâmica do hidróxido Fe(OH),. A precipitação 
deste é tanto maior quanto mais elevado for o pH. 

Como é óbvio, as reacções químicas independentes 
do pH não podem ser representadas num diagrama 
E-pH. 

Representando da maneira descrita acima, num 
mesmo diagrama potencial-pH, as condições de equi- 
líbrio de todas as reacções químicas e electroquímicas 
susceptíveis de se verificarem num determinado sis- 
tema, obter-se-á um emaranhamento de linhas e de 
famílias de linhas que representará de maneira sin- 
tética o conjunto dessas condições de equilíbrio. Tal 
diagrama de equilíbrios electroquímicos permite 
representar as condições de estabilidade ou instabili. 
dade de cada um dos constituintes considerados, e 
isso, em princípio, seja qual for a complexidade do sis- 
tema em estudo. A este propósito veja-se a Fig. 4 onde 
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Fig. k Diagrama do sistema Zn-H0, 


se representa um diagrama muito simples do sistema 
zinco-água considerando apenas as espécies químicas 
Zn, Zn”. Zn(OH), e ZnO;. Como se vê nesta figura, 
verifica-se que o diagrama contém três zonas distintas: 
uma zona de equilíbrio do metal com os seus iões 
Zn'*; uma segunda zona correspondente ao equilíbrio 
do metal com o seu hidróxido Zn(OH), e na qual a 
tensão do metal varia com o pH duma maneira 
independente da natureza do hidróxido; finalmente, 


uma terceira zona de equilíbrio entre o metal e o seu 
complexo solúvel ZnO;. Quer dizer: o zinco com- 
porta-se como um redutor tanto mais enérgico quanto 
mais alto for o pH do meio com o qual está em 
contacto e esta conclusão é válida para a maioria dos 
metais. 

A propósito da influência da água nos equilíbrios 
E-pH, refira-se desde já que nos diagramas de equi- 
líbrio electroquímicos, é vantajoso indicar a provável 
estabilidade dos diversos constituintes em presença 
da água e por conseguinte há que sobrepôr ao dia- 
grama do tipo metal-água o diagrama da água pura 
(a 25ºC e a 1 atmosfera). Na secção que se segue, 
iremos estabelecer e analisar o diagrama da água pura 
após o que ficaremos dotados de ferramenta suficiente 
para poder estudar mais pormenorizadamente um 
diagrama de equilíbrio dum sistema metal-solução. 

Resta-nos a concluir esta secção, fazer uma breve 
referência ao cálculo das constantes Eº e K e colec- 
tâneas onde parte delas se encontram tabeladas. 


Cálculo de constantes-Tabelas 


O uso prático das equações caracterizantes dos 
equilíbrios que se representam nos diagramas E-pH 
implica o conhecimento das constantes Eº e K. 

Quanto ao valor de Eº, as relações: 


Ee = — AGº/nF e Eº=(RT In K)/nF 


mostram que pode ser obtido ou a partir de dados 
energéticos, ou à custa de dados de equilíbrio. 

No primeiro caso, o cálculo de Eº pode fazer-se, 
dadas as propriedades aditivas dos potenciais termo- 
dinâmicos, a partir das energias úteis de formação 
(energias parciais molares) de cada constituinte refe- 
ridos ao mesmo estado-padrão. Porque este estado se 
refere habitualmente a concentrações ou a pressões 
unitárias, aqueles valores confundem-se com os poten- 
ciais químicos dos constituintes na mesma base. 

Assim, a energia global do sistema poderá cal- 
cular-se por uma relação da forma: 


à Gº pode ser calculado ainda por dados termoqui- 
micos uma vez conhecidos 4 Hº e 4 Sº, pois 


AGº=AHº—- TASº 
Igualmente se pode calcular K pela relação: 
AGe=-—-RT Ink 
Em muitos casos os valores de Eº e K encontram-se 
já em tabelas; o mesmo acontece com os potenciais 


químicos padrão. 
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As actuais melhores colectâneas destes valores, 
por mais completas e pela confiança que merece os 
seus autores são: 


— Selected Values of Chemical Thermodynamic Pro- 
perties, Série 3, U.S.A. Bureau of Standards, Was- 
hington, 1954-56. 

— The Oxidation States of the Elements and their 
Potentials in Aqueous Solution, W.M. Latimer, 
2.º ed., Prentice Hall, New York, 1961. 

— Oxidation Potentials, W.M. Latimer, Prentice- 
-Hall, New York, 1952. 

— Mettallurgical Thermochemistry, O. Kubaschewski, 
E. LL. Evans e C.B. Alcock, Pergamon Press, 
London, 1967. 

— Stability Constants, L.G. Sillén, The Chemical 
Society, London, 1964. 

— Contribution to the Data on Theoretical Metal- 
lurgy, XIII High Temperature Heat Content, Heat 
Capacity and Entropy Data for the Elements and 
Inorganic Compounds, ed. K. K. Kelley, U.S. Bur. 
Mines Bull., 584, Washington, 1960. 

— JANAF Thermochemical Tables, Advanced Rese- 
arch Projects Agency Program, U.S. Air Force 
Contract n.º AF 33(616)-6149. 

— Standard Aqueous Electrode Potentials and Tem- 
perature Coefficients at 25ºC, A. de Bethune e 
N. A. Suendenman Loud, Hampel, Skokie, Illinois, 
1964. 

— Atlas d'Équilibres Electrochimiques à 25ºC, M. 
Pourbaix, Gauthier-Villars, Paris, 1963. 


DIAGRAMA DE EQUILÍBRIO ELECTROQUÍ- 
MICO DA ÁGUA 


Sendo a água um dos constituintes principais dos 
sistemas metal-solução aquosa é de toda a conve- 
niência estabelecer e analisar o diagrama de equilíbrio 
electroquímico da água, antes de proceder ao estudo 
dum diagrama do tipo metal-água. 

A presença da água num sistema metal-solução 
aquosa tem consequências muito variadas em virtude 
das suas diversas propriedades físico-químicas que lhe 
permitem actuar, entre outras, das seguintes ma- 
neiras: 


a) Como sistema ácido-base; 

b) Como sistema de oxidação-redução embora 
lento; 

c) Como solvente propriamente dito presidindo às 
dissoluções iónicas e ainda pela sua acção polar 
de hidratação. 


A sua acção mais enérgica é a primeira e deve-se 
ao facto de ser um solvente antiprótico, isto é, tanto 
é dador como aceitante de hidrogeniões. Devido a esta 
dupla acção os iões da água estão sempre presentes 
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num meio aquoso e dadas as propriedades que lhes 
são peculiares, têm acção preponderante sobre os 
estados de equilíbrio das diversas formas sólidas ou 
dissolvidas. 

Quanto à acção óxido-redutora da água pode 
dizer-se que é menos enérgica do que a sua acção 
ácido-base, embora este termo nada represente de 
quantitativo nem seja possível comparar acções tão 
diversas; deverá antes dizer-se talvez, mais restrita. 
Apesar de tudo, a acção óxido-redutora da água 
poderá também influenciar a estabilidade de certos 
constituintes de alguns sistemas, acção da qual resulta 
a decomposição da própria água. 

Convém pois ao estudar um meio aquoso, não 
só atender às influências possíveis da água, como 
fazê-lo de modo a poder avaliar quantitativamente 
essas influências. Daqui a importância do pH das 
soluções, como dos potenciais de oxidação-redução 
da própria água. 

A fim de analisar a acção da água sob o ponto 
de vista ácido-base considere-se então o equilíbrio 
de dissociação existente entre a água e os seus iões 
H* e OH; o qual nos é extremamente familiar: 


H,O = OH + Ht 


A 25ºC, a condição de equilíbrio para esta reacção 
escreve-se: 


Vgts + HoH- = légro 
log (HH) (OH) = log K=—-—-————————— 


1363 
é como His = 0, Hour = — 37595 cal e His = — 56690 
cal [2], resulta log K = — 14,00 ou log (H*) + 
+ log (OH) = — 14,00. 
Visto que pH = — log (H*) tem-se, para 
pH = 7,00 (H*) = (0H) soluções neutras 
pH < 7,00 (H") > (0H) soluções ácidas 
pH > 7,00 (H”) < (0H) soluções alcalinas 


A linha vertical de pH = 7,00 traçada na Fig. 5 separa 
portanto o domínio onde há predomínio de hidro- 
geniões em relação aos iões OH" (soluções ácidas) do 
domínio onde há predominância de iões OH" em rela- 
ção aos 10es H* (soluções alcalinas). 
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É também sabido que a água e seus constituintes 
iónicos pode ser reduzida (com libertação de hidro- 
génio) ou oxidada (caso em que a reacção é acom- 
panhada pela libertação de oxigénio), de acordo com 
as equações: 


2H +2e > H, 
redução em meio ácido (a) 
2H,0+2e >» H,+2 OH 
redução em meio alcalino (a) 
2 HO -» O, + 4H* + de 
oxidação (b) 


Isto conduz às seguintes condições de equilíbrio: 
Reacção 2H* + 2e = H, 


E, = 0,0 — 0,0591 pH — 0,0295 log pg, 


= 0,0 — 0,0591 pH + 0,0295 rH, com rH = 
= — log Px, 


Reacção O, + 4H* + de = 2 HO 


E, = 1,229 — 0,0591 pH + 0,0147 log po, 


= 1,229 — 0,0591 pH — 0,0147 rO, com rO = 
= — log Po, 
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As linhas paralelas a e b de declive —0,0591, traçadas 
na Fig. 6, representam as duas condições de equilíbrio 
para as reacções (a) e (b), com py,=po,=1 atmosfera, 
ou, o que é equivalente, com rH=r0=0, Entre 
estas linhas, as pressões de equilíbrio em hidrogénio 
e em oxigénio são ambas inferiores a 1 atmosfera; 
a região entre estas duas linhas é portanto o domínio 
da estabilidade termodinâmica da água à pressão de 
| atmosfera. Abaixo da linha a correspondente ao 
Pu, = 1 atm. a água sob uma pressão de hidrogénio 
de 1 atm. tenderá a decompôr-se por redução segundo 
a reacção 2H*+ + 2e —» H,. Acima da linha b corres- 
pondente ao po, = 1 atm., a água sob uma pressão 
de oxigénio de 1 atm. tenderá a decompôr-se por 
oxidação segundo a reacção 2H0 => O0,+ 
+4 HH + de. 


Ora, da mesma maneira que a água se pode disso- 
ciar em iões Ht e OH” segundo a reacção HO = 
= OH + H* e é considerada como neutra no ponto 
de vista ácido-base se (H*) = (OH), a água pode 
decompôr-se em hidrogénio e oxigénio gasoso segundo 
a reacção: 


2H0=0,+2H, 


e pode ser considerada como neutra sob o ponto de 
vista da oxidação-redução quando py, = 2po,, quer 
dizer 


log pn, = log Po, + log 2 = log po, + 0,30 
ou se rH = rO — 0,30. 


Se se combina esta relação com as relações E, e E, 

obtém-se, como condição de neutralidade da água sob 
o ponto de vista da oxidação-redução, a 25º C: rH = 
= 27,56 e rO = 27,86. 
Combinando estes valores com a relação pH = 7,00 
que exprime a neutralidade de uma solução aquosa 
sob o ponto de vista ácido-básico, obtêm-se as carac- 
terísticas seguintes para a condição de neutralidade 
absoluta de uma solução aquosa diluída a 25ºC 
(ponto N, Fig. 5): 


pH = 7,00 
IH = 27,56 
rO = 27,86 
E = +0,40V 


Estas características de neutralidade permitem 
dividir um diagrama tensão-pH nas quatro regiões 
representadas na Fig. 5. Esta divisão de M. Pour- 
baix [2] em meios oxidantes ou redutores e ácidos 
ou básicos segundo os valores de pH e do potencial 
de oxidação correspondente, é, sob o ponto de vista 
de coordenação e classificação lógica, muito interes- 
sante; tem também ainda interesse prático no caso 
dos sistemas químicos em que participam oxidantes 
ou redutores enérgicos, e ainda mesmo no caso dos 
sistemas electroquímicos embora as sobretensões pos- 
sam fazer alterar, frequentemente, a ordem por que 
se dão certas descargas. 


EMPREGO E LIMITAÇÕES DOS DIAGRAMAS 
Aplicações em Corrosão 


Os metais com maior interesse económico, tais 
como o ferro, o zinco e o alumínio, pertencem, como 
se pode ver na Série Electroquímica dos Elementos, 
à categoria dos metais «não nobres». A transformação 
desses metais em 10es, seguida pela formação dum 
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sal ou dum óxido, efectua-se quase sempre com liber- 
tação de energia, isto é, trata-se de metais que se 
dissolvem espontaneamente em certas condições favo- 
ráveis, que são usuais na prática corrente (presença 
de humidade, acesso ao oxigénio atmosférico). 

Para obter uma visão de conjunto tão completa 
quanto possivel das reacções de interfase que são 
termodinamicamente exequíveis num caso particular 
determinado, os diagramas de equilíbrios electroqui- 
micos são extremamente úteis. Estes diagramas per- 
mitem prever, como já acima se disse, o que pode 
acontecer a um eléctrodo do tipo metal-solução 
caracterizado por determinados potencial e pH. 

Seja qual fôr o metal, admite-se geralmente que 
um metal é incorroivel quando, mergulhado numa 
solução aquosa inicialmente isenta desse metal, a 
quantidade desse metal que a solução pode dissolver, 
em condições de equilíbrio termodinâmico, não ultra- 
passa 10º átomo-grama/litro (ou seja 0,06 mg/1 para 
o Fe; 0,03 mg/1 para o Ale 0,2 mg/1 para o Pb). 

Por conseguinte, nos diagramas tensão-pH, as 
linhas que correspondem a uma solubilidade de 10* 
átomo-grama/litro (isto é, as linhas de equisolubili- 
dade de cota —6) separam os domínios de corrosão 
e de não corrosão. Neste último domínio, há ainda 
que distinguir entre a imunidade (ou protecção cató- 
dica) que corresponde ao domínio de estabilidade do 
metal, e a passivação que corresponde a uma região 
em que o metal não é estável mas está recoberto por 
um óxido, hidróxido, hidrido ou sal. 

No estado de imunidade o metal é absolutamente 
incorroível porque a reacção de corrosão é energeti- 
camente impossível; por conseguinte, o metal é a 
forma sólida estável. 

No estado de passivação, o metal, embora não se 
encontrando em condições de imunidade, deixa de 
apresentar as propriedades de superfície correspon- 
dentes ao estado metálico solúvel; a forma estável é 
agora um óxido (ou hidróxido, ou sal sólido) inso- 
lúvel. Se este depósito que recobre o metal é aderente 
e não poroso, o metal não sofrerá praticamente 
corrosão dizendo-se que está num estado de passivi- 
dade; se o depósito é poroso pouca protecção confere 
ao metal, isto é, a passivação não implica forçosa- 
mente a ausência de corrosão. 


No caso particular do ferro, por exemplo, o óxido 
que se forma por passivação é, em geral, não poroso 
é protector no caso das soluções isentas de cloretos 
e não protector quando em presença de soluções com 
cloretos; então, no 1.º caso o óxido «passiva» o metal, 
eliminando praticamente a sua corrosão, enquanto 
que no 2.º caso o metal fica «passivado» mas a sua 
corrosão é apenas minimizada. 

Na Fig. 4 estão representadas as condições teó- 
ricas de corrosão, imunidade ou passivação do zinco. 
O estado de imunidade corresponde à zona delimitada 
superiormente pelas rectas representativas dos equilí- 
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brios entre o Zn e os iões Zn”, o Zn e o seu hidró- 
xido Zn(OH), e o Zn e os iões ZnO;. O estado de 
passivação corresponde à zona delimitada lateralmente 
pelas rectas representativas dos equilíbrios entre os 
Zn? e o Zn(OH), e o Zn(OH», e os iões ZnO; , e infe- 
riormente pela recta representativa do equilíbrio 
entre o Zn e o Zn(OH),. Finalmente o domínio 
de corrosão abarca duas regiões: as de estabilidade 
termodinâmica dos iões simples e complexos Zn” e 
ZnO;, respectivamente. Nesta figura a e b repre- 
sentam as rectas limites inferior e superior, respec- 
tivamente, do domínio de estabilidade da água. 
Vê-se, pois, que o zinco é um metal não nobre uma 
vez que as linhas de equilíbrio Zn/Zn” se situam 
fora do domínio de estabilidade da água. Ensaios 
experimentais conduzidos com o zinco ajustam-se à 
interpretação permitida pelo diagrama. 

Duma maneira geral, os metais que resistem à 
água pura, a cerca de 25º C, são aqueles para os quais 
os diagramas potencial-pH correspondentes são tais 
que as rectas pH = 7 só atravessam os seus domínios 
de imunidade ou de passivação aos potenciais com- 
preendidos entre —0,8 e +0,7 V (relativamente ao 
eléctrodo padrão de hidrogénio), potenciais que são 
os únicos possíveis nas condições normais de utiliza- 
ção dos metais de construção, pois que noutras con- 
dições verifica-se a decomposição da água com os 
efeitos consequentes. Estes metais são, para além dos 
seis do grupo da platina, os treze seguintes: berílio, 
alumínio, gálio, índio, estanho, prata, ouro, titânio, 
zircónio, háfnio, nióbio, tântalo e crómio. No caso 
do estanho e do titânio há possibilidade de corrosão 
a potenciais de eléctrodo muito baixos, pelo que con- 
vém evitar acções redutoras e assegurar um bom poli- 
mento da superfície desses metais. 


Aplicação em Electrodeposição Metálica 


Os diagramas potencial-pH permitem-nos obter 
uma ideia sobre as possibilidades de deposição cató- 
dica eficiente de um metal a partir duma solução 
dos seus iões. Para que o problema da libertação 
simultânea do hidrogénio não assuma proporções 
demasiadas (originando depósitos metálicos pouco 
aderentes e conduzindo à precipitação de óxidos ou 
hidróxidos na interfase cátodo-solução) é preciso, em 
princípio, que o domínio de estabilidade do metal 
caia dentro do domínio de estabilidade da água, e 
isto numa gama de pH em que os iões metálicos 
possam existir em concentração elevada na solução. 
Quando tal não acontece, a electrodeposição metá- 
lica só se fará convenientemente se a sobretensão 
associada à libertação do hidrogénio sobre o metal 
for elevada e a sobretensão associada à descarga dos 
iões metálicos sobre o metal for desprezável. Há 
portanto, neste caso, interesse em trabalhar em con- 
dições de electrólise tais que os potenciais reversíveis 


(ou de Nernst) dos sistemas M*/M e H+/HX((M) sejam 
aproximadamente iguais. Repare-se que uma depo- 
sição em meio fortemente ácido tem o inconveniente 
da codescarga do hidrogénio e numa deposição em 
meio ligeiramente ácido corremos o risco de deparar 
com a precipitação de um hidróxido sobre o cátodo. 
A vantagem das soluções alcalinas é pois aparente 
sendo os depósitos metálicos obtidos nestes meios, 
em geral, de muito má qualidade. 

Os diagramas potencial-pH também nos auxiliam 
na questão da deposição de um metal sobre outro 
(metal-base) a partir duma solução que contém iões 
do metal a descarregar. Suponhamos, por exemplo, 
que pretendemos cobrear o ferro. O cobre pode ser 
depositado catodicamente em meio ácido, originando 
um depósito metálico desde que o pH não ultrapasse 
o valor de 2,5 e o potencial catódico seja igual a 
cerca de 0,1V para que não se atinja o potencial 
de descarga do hidrogénio [3]. Acontece que nestes 
meios, o ferro corrói-se a potenciais superiores a 
cerca de —0,7 V, isto é, em banhos ácidos os iões 
Cu” promovem a dissolução do Fe, pelo que são de 
rejeitar. Recorre-se então a banhos cianetados alca- 
linos (de pH 11 a 12) que permitem cobrear o ferro 
devido, não só às características alcalinas desses 
banhos, mas ao facto de que a presença do cianeto 
tem por efeito estender consideravelmente até aos 
potenciais vizinhos de —1V (aos quais o Fe já se 
não corrói), o domínio de corrosão do cobre, por 
causa da formação de complexos cuprosos. Aprovei- 
ta-se para informar que os aniões, como os cianetiões, 
que formam complexos solúveis tendem a aumentar 
as zonas de corrosão, enquanto que os aniões, como 
os sulfatiões, que formam compostos insolúveis ten- 
dem a ampliar a zona de passivação. 


Tal como quando falámos do emprego dos diagra- 
mas E-pH em corrosão, também agora se vê, a pro- 
pósito do seu uso na electrodeposição dos metais, 
que é conveniente recorrer ao diagrama da água pura 
para melhor compreender o alcance das afirmações 
feitas acima. Fica pois mais que justificada a necessi- 
dade de sobrepôr a diagramas de sistemas M-H,O o 
da água pura. 


Limitações dos diagramas 


Os diagramas de equilíbrios electroquímicos dos 
metais constituem uma forma eficiente de represen- 
tação de dados termodinâmicos característicos dos 
sistemas de corrosão. Contudo, apresentam certas 
limitações que devemos ter presente sempre que a 
eles recorramos; assim: 


l. Dizem apenas respeito ao comportamento de 
metais puros (por exemplo, ferro puro), na 
presença de soluções isentas de substâncias 
capazes de formar complexos solúveis (ciane- 


tos) ou sais insolúveis (fosfatos) com esses 
metais. Não se aplicam portanto às ligas me- 
tálicas. 

2. Os diagramas teóricos de «corrosão-imunidade- 
-passivação» que se podem deduzir dos dia- 
gramas de equilíbrios electroquímicos, apoia- 
dos em certas hipóteses de base, aplicam-se 
apenas aos sistemas e condições para os quais 
se constuíram esses diagramas de equilíbrios. 

3. Pressupõem que todas as reacções repesen- 
tadas no diagrama são reversíveis e rápidas. 
Ora, na prática, as reacções entre o metal e os 
seus iões ou entre estes e os óxidos ou hidró- 
xidos do metal, por exemplo, são, em geral, 
irreversíveis, encontrando-se afastadas do seu 
estado de equilíbrio. 


4. Requerem o conhecimento das actividades 
individuais mas como, na prática, estas são 
desconhecidas, substituem-nas nos cálculos por 
concentrações. 

5. Não nos informam acerca da cinética dos 
processos (velocidade de corrosão), visto que 
são diagramas termodinâmicos, ideais. 

6. Aplicam o termo passivação às zonas de esta- 
bilidade dos óxidos e hidróxidos, independen- 
temente das suas propriedades protectoras. 

7. Identificam o pH à superfície do metal com 
o pH no seio da solução. Acontece que o pH 
que se costuma medir (pH do conjunto ou 
pH da solução) é normalmente diferente do 
pH na vizinhança do metal (pH da fracção da 
solução que está em contacto directo com a 
superfície — ou com a fracção da superfície — 
do metal examinado) e, além disso, há também 
uma variação local do pH na superfície metá- 
lica sendo este mais baixo nas áreas anódicas 
(zonas onde se forma um óxido ou hidróxido, 
como é o caso da ferrugem) e mais alto nas 
áreas catódicas (zonas onde o oxigénio dissol- 
vido se reduz, por exemplo). 


A propósito da limitação acabada de apontar, é 
interessante referir que as variações locais do pH, as 
quais se manifestam essencialmente aquando da ocor- 
rência de corrosões localizadas, motivaram uma apa- 
rente discrepância entre os fenómenos observados e 
os previstos pelos diagramas E-pH, o que dificultou 
os estudos tendentes a estabelecer os mecanismos da 
corrosão tipo pitting e de certas formas de corrosão 
sob tensão. Efecivamente a medição de pH junto aos 
metais onde esses fenómenos ocorriam conduziu a 
valores de pH superiores a 8, isto é, de acordo com 
os diagramas correspondentes estar-seiia em zonas 
de passivação sendo difícil explicar a dissolução que 
necessariamente teria de ocorrer nos pits e outros 
locais de dissolução à superfície do metal. A diver- 
gência explicou-se só recentemente [4], quando foi 
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possível medir os pH dentro dos pits e outras caver- 
nas (ou fendas) e observar que, de facto, podem 
verificar-se variações locais da ordem de 5 unidades 
de pH (Fig. 7); assim, Brown e outros [5] determi- 
naram o pH de uma solução exsudada de uma fissura 
propagando-se num aço imerso numa solução neutra 
de NaCl e verificaram que o pH junto ao metal era 
de 6,5 e no interior da fissura foi decrescendo, com 
a sua propagação, para cerca de 3. 


DIAGRAMAS DE SISTEMAS DO TIPO METAL- 
“AGUA 


Baseando-nos em tudo o que atrás dissemos, utili- 
zando os dados quantitativos mencionados, contentar- 
-nos-emos agora com o estabelecimento e a análise 
do diagrama do sistema ferro-solução aquosa, em que 
o pH varia graças à presença dum ácido cujo anião 
se supõe não exercer qualquer efeito complexante 
sobre os iões em solução, nem participar na forma- 
ção de nenhum dos constituintes sólidos postos em 
Jogo. 


Estabelecimento do diagrama 


Para estabelecer um diagrama E-pH é útil pro- 
ceder como se segue: 


a) Estabelece-se um quadro com as fórmulas das 
espécies químicas que se propõe considerar e os valo- 
res dos potenciais químicos padrão uº de todas estas 
espécies. 


No caso presente, e para não alongar esta expo- 
sição consideram-se apenas os seguintes constituintes 
principais: 


Fórmula da espécie 


química nº (cal/mole) 


Formas sólidas: 


Fe O 
Fe;O,; — 177100 
Fe,O, — 242400 
Formas dissolvidas: 
Fe? — 20300 
Fe* — 2530 
HFeoO, — 90627 


b) Escrevem-se as equações das diferentes reacções 
nas quais estes corpos podem tomar parte dois a dois 
fazendo-se intervir eventualmente a água, os seus 
iões e a carga eléctrica livre e; seguidamente cal- 
culam-se as condições de equilíbrio destas diferentes 
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pH mo dedo 809 W 


pM no pi Zoé 


Fig. 7 Arranjo para a medição do ph no interior de um pit 
artificial. 


reacções. No nosso caso, e continuando-nos a limitar 
às reacções principais, teremos: 


1. 2 Fe+t+ 3 HO = Fe,O, + 6H 

2. Fe0O,+ 6H +22 =2 Fe*+3 HO 
3. FeO,+ 8Ht+ 2e=3 Fe* + 4 HO 
4. Fet + 2e = Fe 

5. 3 Fe,O, + 2H* + 2e = 2Fe,O, + HO 
6. Fe;O, + 8H* + 8e = 3Fe + 4H,0 

7. Fe,O, + 2H,0 + 2e = 3HFeo, + Hº 


8. HFeO, + 3H+ + 2e = Fe + 2H O 
9, Fet+ e = Fe” 


log (Fe*) = —0,72 — 3pH 

E = 0,728 — 0,1773 pH — 0,0591 log (Fe*”) 

E = 0,980 — 0,236 pH — 0,0886 log (Fe”) 

E = —0,440 + 0,0295 log (Fe?) 

E = 0,221 — 0,0591 pH 

E = —0,085 — 0,0591 pH 

E = —1,819 + 0,0295 pH — 0,0886 log (HFeO)) 
E = 0,493 — 0,0886 pH + 0,0295 log (HFe0O)) 


e 


E = 0,771 + 0,059] log (Fe). 


c) Traçam-se sobre o mesmo diagrama as dife- 
rentes linhas e famílias de linhas de que se falou 
acima obtendo-se um diagrama de conjunto que 
representa as condições de equilíbrio de todas as 
reacções consideradas. No caso do sistema Fe—H,0, 
considerando como constituintes precipitados o Fe,O; 
e o Fe;O, obtém-se o diagrama representado na 
Fig. 8. Neste diagrama, as linhas a e b dizem res- 
peito, como já se viu aquando do estudo do diagrama 
da água pura, à libertação de hidrogénio e de oxi- 
génio, respectivamente: 


a 2H + 2e=H, 
b. O, + 4H* + 4e = 2H,0 


E = — 0,0591 pH (pu, = 1 atm) 
E = 1,229 — 0,0591 pH (p,, = 1 atm) 


A Fig. 8 representa portanto as circunstâncias de 
equilíbrio termodinâmico do sistema ferro-água, a 
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Fig. 8 Diagrama de estabilidade do sistema ferro - água a 25ºC. 


25º €C, na ausência de espécies complexantes e de 
espécies que formem sais insolúveis com o ferro. 


Interpretação do diagrama 


O diagrama Fe—H,O mostra nitidamente que, 
segundo as condições de pH e de potencial de eléc- 
trodo realizadas, a oxidação do ferro pode conduzir 
a produtos solúveis (lões ferrosos verdes Fe”, 10es 
férricos amarelos Fe”, iões hipoferrosos verdes HFeO”) 
ou produtos insolúveis (magnetite negra Fe;O,, óxido 
férrico castanho Fe,O;,, o qual pode ser diversamente 
hidratado e é o constituinte principal da ferrugem). 
Uma vez que o ferro se pode corroer em presença 
de uma solução quando a quantidade de ferro que 
pode dissolver esta solução é superior a 10% átomo- 
-grama/litro e, por outro lado, que o ferro pode ser 
passivado se se puder revestir com um óxido inso- 
lúvel (por ex. Fe,O,), as linhas que, traçadas na 
Fig. 8, têm cota —6, dividem o diagrama em zonas 
de imunidade, de passivação e de corrosão, tal como 
se representa na Fig. 9, que foi deduzida da Fig. 8 
para o caso em que a passivação é devida à presença 
da película de Fe,O.. 


!- Protecçoo aródica 


2- Protecodo por 
meolinizoços 


3- Protecção cotodiea 


Fig. 9 Corrosão, imunidade e pas- 
sivação do ferro, admitindo 
uma passivação por um fil- 
me de Fed, 


Corrosão 


Como mostra o diagrama, desde que se criem na 
solução em contacto com o ferro condições para as 
quais o Fe não seja estável ele dissolver-se-á ou trans- 
formar-se-á dando-se o que se chama a corrosão. 

Ora verifica-se que o potencial do sistema redox 
2H* + 2e = H, da água é sempre maior que os po- 
tenciais correspondentes à estabilidade do Fe donde 
se conclui que sempre que o ferro é mergulhado 
numa solução não contendo oxidantes enérgicos 
ele pode ser corroido, qualquer que seja o pH dessa 
solução. A maneira como se efectua essa corrosão 
é que é diferente conforme o pH da solução. Assim: 


Soluções ácidas (pH <5)-— Nestas solu- 
ções a reacção directa dos 10es H* da água 
com o ferro é termodinamicamente possível, 
uma vez que o potencial reversível do sistema 
H+/H, é superior ao do sistema Fe?*/Fe. Por 
conseguinte, a corrosão dá-se por dissolução 
do ferro com passagem ao jão Fe” (Fe + 
+ 2H* — Fe” + H,; ponto 1, Fig. 8) e poste- 
rior possibilidade de precipitação de Fe;O, 
para as soluções de maior pH (3Fe* + 
+4H,0 — Fe;O, + 8H* + 2e; ponto 2, Fig. 8). 

Soluções fortemente básicas (pH 5 13,5) — 
Nestas soluções, uma vez que o ferro é sus- 
ceptível de ser complexado pelos iões da água, 
o pH intervém directamente sobre a tensão 
do metal-solução, a qual decresce com o au- 
mento do pH. Por conseguinte, a corrosão 
dá-se por dissolução do ferro e formação do 
ião HFeO) que é estável na água a esses valo- 
res de pH (Fe + 2H,0 — HFeoO; + 3H+ + 20). 

Soluções de pH intermédio — Nestas solu- 
ções a corrosão dá-se com precipitação de 
Fe;O, que rapidamente passa a Fe,O; (passi- 
vação). Há pois possibilidade de reacção di- 
recta entre os H* da solução e o óxido que 
recobre o metal (Fe;O, + 8H' + 8e = 3Fe + 
+ 4H,0). A equação de Nernst correspondente 
a este equilíbrio, admitindo que ele é reversí- 
vel, escreve-se: E = 0,085 — 0,059 pH, isto é, 
o ferro recoberto pelo seu óxido Fe;O, (ou 
Fe,O;, Fe (OH), ou Fe (OH);) comporta-se de 
modo análogo a um eléctrodo de hidrogénio: 
M/MO/H* ou M/MO/0H.. 


Por conseguinte, vê-se que a influência do pH 
sobre a tensão dum eléctrodo resulta da possibilidade 
de se darem os três tipos de reacções de corrosão 
acabados de expôr. No caso do ferro essa influência 
pode pois indicar-se pelas linhas 4'68 sobre o dia- 
grama da Fig. 8. 
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Imunidade 


A imunidade do ferro metálico será pois conse- 
guida sempre que o potencial do ânodo constituído 
por esse metal for obrigado a manter-se, por pro 
tecção catódica, inferior a —(,62 V para soluções de 
pH < 8,9; inferior a valores de potencial variando 
lineamente desde —0,62 a —0,90V quando 8,9 < 
<pH< 13,6; e inferior a valores de potencial va- 
riando lineamente desde —0,90V e segundo a linha 
8' (Fig. 8) para soluções muito alcalinas. 


Passivação 


A dissolução do ferro pode terminar por um efeito 
de passivação, como se viu. Esta pode ter várias 


causas, mas quando se submete o ferro à dissolução 
anódica, a passivação resulta da formação de uma 
película fina, que na região de pH onde a passivação 
é mais fácil (pH — 10 a 12), consiste de Fe;O,. 

No entanto esta passivação é possível num elec- 
trólito de menor pH desde que se deixe o sistema 
evoluir por si. E isto porque a formação de cada vez 
maior quantidade de 16es ferrosos (a que se segue a 
formação de iões férricos Fe”) na região anódica, 
à medida que a densidade de corrente sobe, leva a 
uma alcalinização do meio devido à libertação de H; 
no cátodo, que se dá de preferência à descarga do 
Fe. Criam-se assim condições propícias à formação 
de películas de passivação, particularmente protecto- 
ras se formadas pelo óxido Fe,O, y. 

Também numa electrólise em meio muito alca- 
lino se pode observar a passivação do metal. De facto 
nesta região de pH a dissolução do Fe segundo a 
equação: 


Fe +2HO0 — HFe0O; + 3H* + 2e 


pode levar a uma acidificação do meio sobretudo se 
o anólito estiver convenientemente isolado do cató- 
lito. Esta acidificação, como é evidente a partir do 
gráfico da Fig. 8, leva à formação da magnetite 
Fe;O,. De qualquer maneira, mesmo sem esse iso- 
lamento, a acidificação do meio verifica-se sempre, 
visto que à formação de três hidrogeniões no ânodo 
corresponde a descarga de apenas dois, no cátodo. 


Verificações experimentais 


A fim de verificar e de precisar as conclusões 
deduzidas do diagrama do Fe —H,O, Pourbaix [3] 
fez várias experiências, que confirmaram as previsões 
teóricas. Assim, introduzindo fios de aço em provetas 
com água ou com diferentes soluções aquosas, ele 
observou que: — o ferro se corrói na água destilada 
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(pH = 81; E= —0,486V) com formação de ferru- 
gem castanha; a água da cidade de Bruxelas corrói 
o ferro (pH = 7,0; E= —0,450V); a corrosão é 
aumentada na presença de ácido sulfúrico (1 g/1) 
havendo um desenvolvimento de hidrogénio (pH=2,3; 
E = —0,351 V); a corrosão é diminuida por adição 
de um redutor tal como o bisulfito de sódio (1 g/1), 
sendo nesse caso o metal recoberto de uma camada 
negra de magnetite (pH = 6,4; E= —0,372V); a 
corrosão pode ser suprimida por adição de cromato 
de potássio (pH = 8,5; E = +0,235 V). Noutra série 
de experiências, Pourbaix introduziu fios de aço, 
ligados por ligação simples a fios de outros metais, 
em provetas contendo água de distribuição urbana, 
e verificou que o par com cobre ou com platina 
provoca um acréscimo da velocidade de corrosão e 
o par com zinco ou com magnésio suprime geral- 
mente a corrosão: 


Par Fe-Cu; pH = 7,5 E=-0,445V —> 
O ferro corrói-se 

Par Fe-Pt: pH = 7,5 E=-0444V —s 
O ferro corrói-se 

Par Fe-Me: pH = 7,5 E=-090V —— 


O ferro não se corrói 
Par Fe-Zn: pH = 7,5 
O ferro não se corrói 


E = -—0,6900V —> 


Assinalemos a finalizar que, de acordo com as 
considerações feitas a propósito do diagrama do ferro, 
é possível protegê-lo contra a corrosão por três pro- 
cessos, como se indica na Fig. 9: —ou se desloca o 
ponto representativo do metal para baixo (protecção 
catódica), para cima (protecção anódica ou passivação), 
ou para a direita (protecção por alcalinização). A título 
exemplificativo, introduzamos uma vareta de ferro nu- 
ma solução ácida (pH = 3); nesta solução o ferro, em 
condições de abandono, corrói-se, assumindo uma 
tensão mista mais alta que a sua tensão reversível 
(ponto 1, Fig. 8)*. Suponhamos agora que preten- 
demos proteger catodicamente o ferro, podendo para 
isso fixar o seu potencial num valor correspondente 


* Repare-se que pela mesma razão que o sistema redox 
2H+ + 2e = H, é capaz de elevar a tensão do metal para o 
ponto 1, também o sistema redox O, + 4H+ + de = 2H,0 a 
poderá elevar para um ponto (ponto 3, Fig. 8) algures abaixo 
da recta limite superior do domínio de estabilidade da água, 
É por isso que no caso dos sistemas do tipo metal nobre- 
solução aquosa, em que as linhas de equilíbrio do metal com 
o seu catião se situam dentro do domínio de estabilidade 
da água, é conveniente desarejar a solução se se pretende 
fazer coincidir a tensão efectiva com a que traduz o equilíbrio 
reversível do metal com os seus iões. O equilíbrio 0,/H,0 
não é no entanto relevante, especialmente no caso dos metais 
comuns em meios ácidos, em virtude da fraca solubilidade 
do O, na solução e porque a corrente de permuta para o 
equilíbrio 0,/H,0 é muito menor que a corrente de permuta 
para o equilíbrio H+/H, (no caso do Fe é 10º vezes menor). 


É 5 | R Mm ph 
Fig. 10 Diagrama potencial-pH do chumbo 
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Soluções salinas E 
Suítotos, corbonotos, bscorá Fell. mtrotos, comutos 
NoCiO, puro ds NOCIO, com Noll 
Cl, seco ou húmido,< DOC El > 100€ 
Br, seco, temp. ombiente Br húmido ou o oito 
502. SD. Ma5 temp. HE 


a um dos pontos 4, 5 ou 6 da Fig. 8. É evidente 
que o potencial a utilizar deveria ser o correspondente 
ao ponto 4; acontece que na prática é difícil definir 
com precisão o ponto 4 (coincidente com o potencial 
de equilíbrio reversível entre o metal e os seus iões), 
pelo que se deve escolher um ponto ligeiramente 
abaixo desse (ponto 5). O ponto 6 não deve ser 
usado, pelas seguintes razões: —a protecção cató- 
dica fica mais cara pois há um maior consumo de 
energia (note-se que o AE=E— E, é muito 
maior que o AE=E,-—E,); a potenciais tão 
baixos como o E, há franca libertação de hidro- 
génio a qual pode conduzir à fragilização ou empola- 
mento pelo H;, do ferro com o consequente empo- 
brecimento das suas propriedades mecânicas; o desen- 
volvimento de hidrogénio pode levar a uma enorme 
alcalinização da solução que contacta o ferro, evo- 
luindo o sistema do ponto 6 para um ponto 7, isto é, 
o ferro deixa de estar no seu domínio de imunidade 
para cair no seu domínio de corrosão. No caso de 
se fixar o potencial num valor correspondente ao 
ponto 5, a libertação de hidrogénio também pode 
deslocar o pH local, mas mesmo que isso aconteça 
o sistema cairá no seu domínio de passivação (ponto 
8, Fig. 8) que, em princípio, assegurará uma protecção 
perfeita ao metal. Como se vê, a protecção catódica 
não é isenta de perigos, o mesmo se podendo dizer 
das outras vias de protecção preconizadas. Este 
assunto, no entanto, sai do âmbito desta exposição. 


CONCLUSÕES 


Os diagramas de equilíbrios electroquímicos devem 
considerar-se principalmente como esquemas que 
indicam, para diferentes condições de potencial de 
eléctrodo e de pH, as reacções que são teoricamente 
possíveis e as reacções que são teoricamente impos- 
síveis. Uma vez que se trata de diagramas termodi- 
nâmicos há que usá-los com plena consciência do seu 
significado e suas condições de validade. Estes diagra- 
mas exprimem essencialmente o que é energetica- 
mente possível ou energeticamente impossível; ora, 
numerosas reacções electroquímicas são irreversíveis, 
isto é, não se produzem efectivamente (ou só se 
verificam muito lentamente), embora sejam energeti- 
camente possíveis. Por conseguinte, estes diagramas 
não indicam necessariamente o que se produz de 
facto; uma resposta formal e definitiva a este pro- 
pósito implica estudos experimentais cinéticos, os 
quais podem consistir, por exemplo, na determinação 
de curvas de polarização intensiostáticas ou poten- 
ciostáticas. 

Apesar das suas limitações, os diagramas são muito 
úteis no estudo de variadíssimos sistemas encarados 
sob os aspectos mais diversos. Aliás, a sua utilidade 
recrudesce quando se modificam tendo em atenção o 
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Fig. Diagrama potencial /pH do ferro. 
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Fig. 12 Diagrama potencial-pH do níquel. 


Resistência à corrosão do níquel em diversos meios 
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Fig. 1E Diagrama potencial-pH do cobre. 


Resistência à corrosão do cobre em diversos meios 
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Fig.13 Diagrama potencial-pH do crômio. 


Resistência à corrosão do crômio em diversos meios 
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Fig. 15 Diagrama potencial-pH do zinco. 
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comportamento real dos metais em meios ácidos, 
alcalinos, oxidantes e redutores. Obtêm-se então dia- 
gramas empíricos que fornecem uma ideia mais 
verdadeira do comportamento corrosivo dos metais e 
que têm, por conseguinte, maior interesse prático, 
Nas Figs. 10 a 15 representam-se os diagramas empí- 
ricos E-pH para os seguintes metais técnicos: Pb, Fe, 
Ni, Cr, Cu, Zn, em meio aquoso não complexante, 
a 25º C. Ao lado de cada diagrama também se indi- 
cam as características corrosivas do metal correspon- 
dente em diversos meios electrolíticos. 

Resta observar que os diagramas de Pourbaix que 
foram traçados para a temperatura de 25º C, passa- 
ram a ser publicados para temperaturas mais elevadas, 
para o que muito contribuíram as expressões teóricas 
estabelecidas por Bethune e outros [6, 7] no princípio 
da década passada. A aplicação do método de Pour- 
baix ao estabelecimento de diagramas semelhantes 
para sistemas metal-sal fundido começou a fazer-se 
há cerca de 20 anos, reportando-se o primeiro trata- 
mento termodinâmico, publicado por Edeleanu e 
Littlewood [8], ao estabelecimento dos diagramas de 
equilíbrio para os sistemas do tipo metal-mistura de 
cloretos fundidos, a 600º C. Em 1966, Perichon [9] 
iniciou a publicação de diagramas de equilíbrios elec- 
troquímicos para sistemas do tipo metal-solvente 
orgânico. Mais recentemente, Pourbaix [10] construiu 
diagramas análogos na presença de meios gasosos. 
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SIMULAÇÃO DE UMA OPERAÇÃO UNITÁRIA 
DA INDUSTRIA QUÍMICA: 
EVAPORAÇÃO DE MULTIPLO EFEITO 


JOÃO FERNANDO PEREIRA GOMES (*) 
CARLOS MANUEL ALVAREZ PEREZ (**) 


RESUMO: 


Apresenta-se um programa de cálculo automático 
em linguagem FORTRAN-IV, para a determinação 
do número óptimo de efeitos de uma unidade de 
evaporação de múltiplo efeito. 

O programa elaborado é geral, e a sua validade 
só é restringida pela validade dos dados fornecidos. 

O programa foi testado com vários exemplos, 
sendo particularmente adequado à desalinização; e 
encontra-se disponível no Centro de Informática do 
Instituto Superior Técnico. 


ABSTRACT: 


This paper presents a computer program using 
FORTRAN-IV language, that determines the opti- 
mum number of effects of a multiple-effect eva- 
poration system. 

The vality of the program is restricted only by 
the vality of the data supplied. 

The program was worked out with various exam- 
ples and applies to desalination. It is available at 
Centro de “Informática do Instituto Superior Téc- 
nico. 


(*) Eng.” Químico-Monitor/Investigador do Grupo de 
Estudos de Processos de Separação-I.S.T. e Assistente Esta- 
giário do L.N.E.T.l. — Dep. de Energias Convencionais. 

(**) Eng.” Químico-Assistente de Produção da HO. 
VIONE, Sociedade Química, Lda. 
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| — PROBLEMA PROPOSTO: 


Pretende-se determinar o número óptimo de efei- 
tos de uma unidade de evaporação de múltiplo efeito 
como o representado na figura 1. É esta a situação 
que mais ocorre na indústria química, constando 
assim de alimentação apenas no primeiro efeito e 
para a frente, havendo uma descarga de cristais em 
cada efeito. 

Por simplicidade de cálculo iremos, no entanto, 
admitir as seguintes aproximações: 


i) Considera-se que a temperatura e a pressão 
das correntes de saída de cada efeito, são as 
mesmas do que as verificadas no corpo do 
evaporador. 


1) Consideram-se desprezáveis as trocas de ener- 
gia com o exterior, nomeadamente as perdas 
por radiação. 


iii) Consideram-se desprezáveis os calores de cris- 
talização. 


iv) Considera-se desprezável a pressão hidrostá- 
tica, e a elevação do ponto de ebulição das 
soluções. 


v) Considera-se que o vapor alimentado a cada 
efeito é vapor saturado, que é condensado sob 
a forma de líquido saturado, após se ter efec- 
tuado a transferência de calor, saindo assim 
à temperatura de saturação. Com isto, des- 
prezam-se sobre aquecimentos e sobrearrefeci- 
mentos, e supõe-se que o calor cedido pelo 


Vi, T; f P, 


A 


ae o E o rm 


V 
Via a 
Ti Tamos 
fia Pio 
q 
CONBÍNVEA PO com nando aa 1 COM BÉMEA Do r 
ALinEa ração e 
CMisrais ; CRISTAIS a 
CRIETAIS ! 
E c. Cm 
= 
Fig. 1 — Sistema considerado 
fluido à solução é apenas o calor latente de 2.1.3 — Balanço global ao efeito À: 
vaporização. 
po da L=L+HV,+C, 
vi) Considera-se que o calor específico do licor in i i i 
não varia significativamente no intervalo de V PT 
temperaturas considerado. ag A tá dh 
| uno Vs 
A validade dos cálculos não vai vir afectada por p 
estas aproximações que são as habituais em casos ; ; 
4 


industriais. 


2 — RESOLUÇÃO DO PROBLEMA: 


O problema é resolvido com recurso a balanços 
mássicos e entálpicos a um efeito genérico i (que se 
representa esquematicamente na figura 2), e ainda 
utilizando as equações de transferência de calor, 
como se verá. 


2.1 — Balanços mássicos 


2.1.1 — Balanço global ao soluto: 
n-1 


Fx =x € + xl 


2.2.1 — Balanço ao soluto no efeito à: 
L,, AL; — L; XL; + C; do; 


Como x, está relacionado com a água de cristaliza- 
ção retida pelos cristais, e é geralmente um dado do 
problema, podemos tomar x. = x. = constante, pelo 
que se obtém: 


LicaR 


Li=4 
Th 4 cfieru 
x 
4-1 Ci 
Ti 
Fig. 2 — Efeito genérico i Ec; 


2.2 — Balanços entálpicos: 
O balanço entálpico ao efeito i escreve-se: 
La Bis TP VatS Lihi + VT 


ED Ms h... “e C, Dcrist | 
sendo h e H entalpias específicas como segue: 
entalpias do licor = h, = Cp T; 
entapias dos cristais = hei; = C'pT, 
entalpia do vapor = H, = A + CpT, 


entalpia do licor = h. | = CpT,, 
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e sendo o estado de referência o das tabelas de vapor. 


2.3 — Balanço combinado de entalpia e massa: 


Combinando a equação anterior do balanço ener- 
gético com os balanços de massa já referidos, obtém- 
-se a seguinte equação: 


XL 


Liz Ai (1 — ) + Lim - dh — A — Cp AT; + 


Xi 


AL. 
a a (Mt A; + T; (Cp — C'p)) | + L; 
Xe 


X, 
| — MT, (cp — cp) |=5 
4 


sendo AT, = T, — Tu. 


Esta é pois a equação do balanço combinado de massa 
e entalpia para um efeito genérico i, exceptuando o 
primeiro efeito, que tem que ser resolvido à parte, 
dado que é alimentado pelo vapor da rede e por 
alimentação fresca. 


Neste caso tem-se: 


F he + Vo E se L; h,, + V; H, + CC, Ncrist; 


em que: 


entalpia específica da alimentação = hp = CpT, 
calor latente cedido pelo vapor da rede = H, = A, 


e as outras grandezas são as já referidas anterior- 

mente, sendo o índice i feito igual a 1. 
Introduzindo todas as relações de balanço mássico, 

desenvolvendo e simplificando, obtém-se então: 


2“F 
Vo=—— | N+Cp (T—TD+ 


Xr 
ha X 


AHTUG-C) | 


E 


at E XL; | 
à |—h+ Es E FT =) 


Repare-se que, com as equações anteriores, se tem 
um sistema de equações lineares em que as incógnitas 
são L,, isto é, o caudal do licor que deixa cada efeito. 
Este sistema de equações pode ser escrito sob a 
forma matricial sendo assim mais simples a sua 
resolução. 
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2.4 — Equações de transferência de calor: 


Além das equações até aqui apresentadas é ainda 
possível escrever as equações de transferência de 
calor que, para o efeito i, se traduz em: 


q = Vi My = UA; (Tn — 1) 


2.5 — Metodologia de resolução: 


A solução do sistema de equações exige que a 
distribuição de temperaturas, ao longo do sistema, 
seja conhecida. Como ponto de partida, admite-se 
uma queda uniforme de temperaturas em cada efeito. 
Assim sendo, a queda global de temperaturas, isto 
é, a diferença entre a temperatura do vapor e a tem- 
peratura no último efeito, é dividida igualmente pelo 
número de efeitos, como segue: 


AT, = A Toxbai/N 


Conhecida a temperatura em cada efeito, proce- 
de-se do seguinte modo: 


1) Determinam-se os caudais do licor que deixa 
cada efeito, L;, resolvendo simultaneamente as 
equações apresentadas. 


2) Calcula-se o vapor alimentado ao primeiro 
efeito. 


3) Usa-se a equação de transferência de calor 
para determinar a área em cada efeito. 


4) Se as áreas não são iguais, a menos de um 
critério definido (por exemplo 1 metro qua- 
drado). as temperaturas são recalculadas, pela 


expressão: 
e Pacioe 
a Lavi di À pterdor E Ee | 
média 
e faz-se: 


N-1 
AT, = A Tatobal e AT, 
i- no 


vo 


5) Repetem-se os passos 1, 2 e 3, até se obterem 
áreas Iguais. 


6) Determina-se o custo da unidade por critérios 
de optimização e de rentabilidade. 


7) O cálculo pode ser repetido para um número 
qualquer de efeitos. 


3 — PROGRAMA DE COMPUTAÇÃO: 


Elaborou-se um programa de computação, em 
linguagem FORTRAN-IV, para a resolução do pro- 
blema proposto, de acordo com a metodologia atrás 
indicada. 

No APÊNDICE I apresenta-se o fluxograma deste 
programa. 

O programa encontra-se gravado em disco do 
Centro de Informática do IST e está à disposição 
dos utilizadores, principalmente dos alunos do curso 
de Engenharia Química. 

Na biblioteca do Departamento de Engenharia 
Química do IST, encontra-se ainda um manual de 
utilização deste programa, assim como uma listagem 
e uma descrição pormenorizada do mesmo. 

O programa foi testado com êxito num exemplo 
concreto de desalinização, utilizando uma unidade de 
evaporação de 10 efeitos, que se encontra também 
descrito no já referido manual. 

Apresenta-se, ainda em apêndice, parte do «out- 
put» obtido nesta aplicação. 

Para este programa, a memória necessária é de 
54 K Bytes, existem cerca de 300 instruções executá- 
veis, e o tempo de passagem típico, em computador 
IBM, para o exemplo referido é de 3 segundos e 20 
décimos. 


NOMENCLATURA: 


A, = área de transferência de calor no efeito i 

- massa de cristais retirados no efeito i 

C, = calor específico a pressão constante, do licor 

C' = calor específico a pressão constante, dos cristais 
F = caudal de alimentação 

h, = entalpia específica do licor que deixa o efeito i 
h,,= entalpia específica do licor que deixa o efeito i 


Roristi 

H, = entalpia específica do vapor que deixa o efeito i 

he = entalpia específica da alimentação 

L, = caudal de licor que deixa o efeito i 

N = número de efeitos 

P, = pressão no efeito i 

q; = quantidade de calor transferida no efeito i 

T, = temperatura no efeito i 

T; = Temperatura da alimentação 

U, = coeficiente global da transferência de calor no 
efeito i 

V, = caudal de vapor que deixa o efeito i 

V, = caudal de vapor da rede alimentada ao primeiro 
efeito 


O 
II 


= entalpia específica dos cristais que deixam o 


APENDICE I 


X.;= concentração de sal nos cristais que deixam o 
efeito i 
X; = concentração de sal na alimentação fresca 
X = a a de sal no licor que deixa o efeito 
; |] 
A T, = queda de temperatura no efeito i 
A Toiobal = queda de temperatura global 
À; = calor latente de vaporização do vapor que deixa o 
efeito i 
À, = calor latente de vaporização de vapor da rede 
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APÊNDICE II - Parte do "output" obtido num caso de desalinização 
para um máximo de 10 efeitos : resultados finais 
para 6 efeitos. 
CALCULOS PARA O SEGUINTE NUMERO DE EFEITOS + & 
TEMPERATURAS EM CAba EFEITO (GRAUS CENTIGRADOSO: 


GE 


EM 


[3 £3T ily.das lig,IJ44 Sa. al O acl UE a 


CONCENTRáCloS DF SAL NO LICOR QUE DEIXA Cal EFEITO <FRACCAO MASSIC 


da DDS A 2B Ba LEI o RE Magra DE ro 


CAuDAL DE LICOR QUE DEIXA LADA EFEITO cko/Ro: 

L/'Bs. e ao 13239: 1103. 141 3?. Sr09. 

COUDAL DE LAPOR QUE DEIZA CADA EFEITO ckG/Ho: 

2h47,385 pedra, Fra 206. 121 1óFB.BZs 197? 1190 «3253 
TEMPERATURA DO VAPOR DA REDE cORAUS CENTIGRADOS:: 17B 

LAUDAL DE VAPOR ALIMENTADO &U FRIMEIRO EFEITO <KG/H)2 v7a4P7.is 
AREA DE TRANSFERENCIA DE CALOR EM CADA EFEITO cMgo: 


Dara ope 107 L% «475 1>. 2HR pd Dr 


ne, 


FRESSAU EM CADA EFEITO CKECMZo: 


tu 
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MODULAÇÃO DE TIPO OQPSK 
COM DIAGRAMA DE OLHO 


PREDETERMINADO 


ANTÓNIO GUSMÃO e NELSON ESTEVES, 
CAPS, IST 


RESUMO 


É apresentada uma técnica de modulação de tipo 
OOPSK («Offset Quaternary Phase Shift Keying») 
em que o sinal a transmitir é sintetizado com base na 
forma do diagrama de olho dos sinais binários modu- 
lantes. 

Mostra-se que este método original de síntese de 
sinais permite englobar num mesmo formalismo 
várias classes de sinais conhecidas e outras inteira 
mente novas e, além disso, oferece uma base de pro- 
jecto flexível e eficiente para condições típicas da 
comunicação digital via rádio. 


ABSTRACT 


An OQPSK-type modulation technique is descri- 
bed where the signal to be transmitted is designed 
on the basis of a prescribed eye diagram . 

Several modulation formats can be described by 
the proposed technique, some of which have not been 
published before and are shown to provide a flexible 
and efficient design basis under conditions frequently 
met in digital radio. 


1. INTRODUÇÃO 
Em comunicação digital via rádio (CDVR) são 


bem conhecidas as vantagens das técnicas de modu- 
lação que dão origem a sinais de envolvente constante 
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ou quase constante e fase contínua. Entre estas téc- 
nicas, têm merecido uma atenção muito especial as 
que, para efeitos de detecção coerente, permitem usar 
um desmodulador com a estrutura simples do desmo- 
dulador convencional dos sinais OQPSK [3-9]. 

Quase todas estas técnicas de tipo OQPSK são 
baseadas no «formato OQAM» («Offset Quadrature 
Amplitude Modulation») ou no «formato CPM» 
(«Continuous Phase Modulation») [1,2]. Estes for- 
matos estão definidos genericamente no Quadro I: 
fo representa a frequência da(s) portadora(s) e q, uma 
fase constante; h é o chamado «índice de modulação»; 
x(t) e g(t) são impulsos na banda de base; os a;'s é by's 
representam símbolos Mr-ários estatisticamente inde- 
pendentes e equiprováveis, o mesmo se podendo dizer 
dos u,'s; finalmente, T é a duração de cada bit à 
entrada do emissor. 

Neste artigo é apresentada uma nova classe de 
sinais de tipo OQPSK caracterizada por quatro fun- 
ções geradoras y;(t), |) = 0,1,2,3. 

Estas funções estão relacionadas de forma simples 
com o diagrama de olho dos sinais binários que mo- 
dulam duas portadoras em quadratura: a Fig. 1 ilustra 
esta relação. 

O método de modulação proposto — que é des- 
crito na Secção 2 — será designado abreviadamente 
por EDB-OQPSK («Eye -Diagram-Based -OQPSK)»). 

Na Secção 3 é estudada a alteração introduzida 
num sinal EDB-OQPSK por um amplificador não 
linear do tipo correntemente utilizado em CDVR. 

Na Secção 4 mostra-se que o formato proposto 
permite descrever várias classes de sinais OQAM e 
CPM conhecidas, e que além disso permite gerar 
sinais inteiramente novos. 
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Na Secção 5 referem-se algumas propriedades úteis 
dos sinais EDB-OQPSK, dando-se especial relevo à 
eficiência espectral. 


2. Estrutura dos sinais EDB-OQPSK 


Na Fig. 2 está representado, muito simplificada- 
mente, um emissor capaz de gerar qualquer sinal 
EDB-OQPSK a partir de um sinal binário de entrada 
que, por conversão série-paralelo, origina o par de 
sinais I(t) e Qu(t). Estes sinais são da forma 


e | t+T-—2KI 
I(t) = e ax TECt (A T ) 
| t+HT-—2KT 
Qt) = A by recr( “92 ) (1) 


com os a;s e b;'s representando símbolos binários 
independentes e tomando os valores +l e —1 com 
igual probabilidade 1/2 (rect(x)=1 se |x|<1/2 e igual 
a zero no caso contrário). 

O modulador propriamente dito é um modula- 
dor convencional de portadoras em quadratura, pelo 


(B) 


Fig. 1 — Alguns exemplos de diagramas de olho (A) 
e funções geradoras correspondentes (B) 
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Y, ltiz ya 1) 


AE 1,02 =1 


QUADRO 1 


— Formato OQAM 
s(t) = io ar x(t—2KT) cos (Qriit+o,) — 
— do be x(t—T—2KT) sen (Qrft+g,) 
ax, bx = +1, +3,..., +(M-— 1) (M par) 
— Formato CPM 


s(t) = A, cos | eta + hr us qgt—KT) + % | 
K 


us = =1,+3,.., <(M-1)(M par) 


t 
q() = f g(r) dr 
— (O) 
+ 00 
f g(Ddt=1 
— 00 


Sinal 
NRZ 
(Ritmo 1/T' 


Fig.2 — Representação simplificada do emissor 
que o sinal EDB-OQPSK pode ser escrito na forma 


s(D=I(t) cos (Qrft+Ho, — QU!) sen (Qrft+Ho) (2) 


sendo f, a frequência da portadora e q, uma fase 
constante arbitrária. 

Os sinais na banda de base, It) e Q.(t), são defi 
nidos do seguinte modo: 


E(t) = Ls Essm(t) 


Qt) = E Qumít) (3) 


com Ism(t) e Qum(t) nulas fora do intervalo [mT, 
(m-+1) T]. Estas funções I;m(t) é Qum(t) dependem 
dos quatro bits apresentados ao codificador durante 
o intervalo [mT, (m+1) T] e também de quatro fun- 
ções geradores y;(t), )=0, 1, 2,3, supostas nulas fora 
do intervalo [0, T], que caracterizam um dado código 
da classe proposta. 


O Quadro II resume o processo de codificação: 


QUADRO II 


CODIFICADOR 


qu to a, q qu ci; a 


EDB- OQPSK 


Se desejarmos um sinal EDB-OQPSK com fase e 
envolvente contínuas, fácil é demonstrar, pela aná- 
lise do processo de codificação, que as funções y;(t) 
também devem ser contínuas, para O<t<T, e que 


yíXD) = y(1)=0 


YoCT-) = yXT”) 
YO) = yolO) 
y(0") = 40º) (4) 


Numa realização possível do codificador EDB- 
-OQPSK — recomendável pela sua simplicidade e 
flexibilidade — a leitura de duas memórias digitais, 
com um conteúdo de acordo com o conjunto (y;(t), 
]J=0,1,2,3) escolhido, permite gerar as sucessivas 
amostras, codificadas em binário, dos sinais It) e 
Qt). Depois da conversão D/A os sinais It) e Qt) 
são aplicados ao modulador, cuja saída é o sinal 
EDB-OQPSK. 

Observe-se que o esquema de codificação indi- 
cado, correspondente ao formato EDB-OQPSK, asse- 
gura que 


Intervalo [20T, (2n+1) T] 


la an(t) 
— Anyát—2nT) 
“any(t—2nT) 
| an yt—2nT) 
a, Y(t—2nT) 


Oss 2n(t) 


ba yo[(Qn+1) T=) 
— bya[Qn+D) T>t] 

b, yi [Qn+I) T—t] 
byslQn+DT=H — 


Intervalo [(Qn+ DT, (Qn+42) T] 


Es, ans(t) Oss ansi(t) 


ba Yolt—(2n+1) T] 
ba yi[t—(Qn+1) T] 
b, y2[t—(Qn+1) T] 
ba ys[t—(2Qn+1) T] 


art YoLQn+D T=] 

As Yo [(Qn+2) T—t] 
ans Yi [Qn+2) T—t] 
any [(2n+2) T>t] 


nl 


s (t] 


Ne; 

Es 

“a 

A=1,0=1 A=/V2a=W2 A,=08,0=08 


Fig. 3 — Exemplo de sinais I,(t) e Q,(t) 


E(CZKkT) = ar yi(0*) ou ar y(0*) 

Q[(2k+1)T] = br yi(0') ou by y(0*) (5) 
uma propriedade importante que sugere a aplicabi- 
lidade, em detecção coerente, de desmoduladores con- 
vencionais «tipo OQPSK», desde que yi(0) e y;(0') 
sejam ambos positivos (ou ambos negativos). 

Na Fig. 3 apresentamos, como exemplo, os sinais 
I(t) e Qt) no intervalo [0, 8T], resultantes das se- 
quências de bits 

a=+l,a=+l,a=-l,a=+l,a=+1 
bi;i=-—1,bo= 1, b=-—l1, b=-—1, b=—1, b= +] 


e do conjunto de funções geradoras 


Yot) = Ag—(Ag—u) sen? ( EE) 


2T 
Yi) = As cos( ZT) 
vt) = 1— (1a) sen? [) 
yi(t) = cos 2T) 
PERDE (6) 
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cujos parâmetros A, e « satisfazem as condições 
O<A<1, O<a< 1. Estas funções e os diagramas 
de olho correspondentes estão representados na Fig. 1, 
para diversos valores de A, e q. Deve notar-se que 
não há uma correspondência biunívoca entre o dia- 
grama de olho e o conjunto ordenado das funções 
geradoras: em cada caso existem duas ordenações 
possíveis satisfazendo (4). 


3. Amplificação não-linear dos sinais EDB-OQPSK 


A envolvente de um sinal EDB-OQPSK pode exi- 
bir flutuações, mas também pode ser rigorosamente 
constante. Por observação do Quadro II podemos 
concluir que a envolvente é constante se, e só se, 


YLD+HyYXT—) = yXD+yXT—U = 
= y(t)+yT—t) = Constante (7) 
para O<t<T. 


Consideremos agora que um sinal EDB-OQPSK 
atravessa um dispositivo não-linear idêntico aos vul 
garmente utilizados em CDVR, para efeitos de am: 
plificação de potência [10]. 

A entrada da não-linearidade o sinal é 


s(t) = It) cos (Qrft+v,))— QU) sen (Qrít+Ho,) = 


= A(t) cos [rft+o(t)+gç] (8) 


com a envolvente A(t) = V Pt) + QD 


ea fase q(t) = tg! [QUU/IC0]. 


O sinal na saída pode ser dado por 


st) = fr [A(O] cos [2nft+o() + +gw [A(O] | 
(9) 


representando fw(.) e gni(.) as funções de «conver- 
são AM/AM» e «conversão AM/PM», respectiva- 
mente. 


É possível provar que o sinal na saída do dispo- 
sitivo não-linear é um novo sinal EDB-OQPSK 
(idêntico ao de entrada se a envolvente deste for 
constante), contanto que gw [A(t)] seja constante, o 
que sucede, nomeadamente, com os limitadores ideais 
(«hard-limiters»). Admitindo que as quatro funções 
características do sinal de entrada são as do conjunto 
(yo), ]J=0,1,2,3), as funções de saída podem ob- 
ter-se das primeiras da forma seguinte: 


Yo (t) 
vo +ydT—o 
x fa ( VEO ET (TO ) 


Ya (t) 


VyX D+ 1) 
X fa ( VEOFETO) 


HO: 


E 


= Ya (1) 
y o (t) — e 
VR FT) 
ya (1) 
so = | 


Vy0+yXT—) 


X fn (vEO FRA ==.) 
(O<t<T). (10) 


4. Algumas subclasses de sinais EDB-OQPSK 
4.1] — Sinais enquadráveis no formato OQAM 


Consideremos um sinal OQAM (ver quadro T) 
com M=2, x(t)=x(—t) e x(t)=0 para [|t>2T. 

Não é difícil provar que um tal sinal pode ser 
encarado como um sinal EDB-OQPSK. As funções 
características são, neste caso, iguais duas a duas: 


Y(D= yít) = x(t) + x4—27T) 


y(D= ylt) = x(t) — x(t—2T) 
(O<t<T). (11) 


Se o impulso x(t) se apresentar restringido na sua 
duração de modo a termos x(t)=0 para |t|>T, as 
quatro funções características serão todas iguais: 


y(t) = x(t), J = 0, Um Es 3 (O<t<T) (12) 


4.2 — Sinais enquadráveis no formato CPM 


Consideremos agora um sinal CPM (ver qua- 
dro I) com h=1/2, M=2, g(D)=g(3T-—t) e g(t)=0 
para té [0, 3T]. Pode provar-se que um tal sinal é 
idêntico a um sinal dado por (2), (3) e pelo Quadro II 
desde ques 
Eh » u,=0 ou x (mod. 27), k par 


é 1i< k 


= + 7/2 (mod. 2x), k ímpar 


b=sen (> | u) (13) 
i< 2k-—1 


=> 


E] 
e q(t) e | g(x) dr se relacione com o conjunto 


— (DO) 


(y(D, 1=0,1,2,3) do modo seguinte: 
Yolt) = Avsen] 5 (s GErZD — ak TD +06) | 


vil) = As sen 5 Gt+2T) — qt+T) — a() )| 


| 


via) E Ansen] 5 (atc+2n + Qt+T) — a(o) 
vio = Às sen) 3 (at+20 + + at+T) + ay )| 


(O<t<T) (14) 


Se o impulso g(t) se apresentar restringido na sua 
duração, de modo a termos g()=0 para tá [T,2T] 
— 0 que equivale a CPM com «resposta total» [1] — 
conclui-se facilmente que as quatro funções gerado- 
ras são todas iguais: 


vit) = A, cos | a q(t+ D |] 511=0,1,2,3(0<t<T) 
| (15) 


significando que um tal sinal CPM é também OQAM. 
É o que sucede, por exemplo, com os sinais MSK [5]: 
neste caso g()=1/T, para T<t<2T, anulando-se fora 
deste intervalo, e portanto 


Tl 
y;i(t) = Açcos = or ). 


4.3 — Outras subclasses 

Uma subclasse de interesse é a que apresenta as 
funções características dadas por (6). Com A, e « 
adequadamente escolhidos, podemos ter sinais com 
envolvente aproximadamente constante e que, por 
tal motivo, o atravessamento de não-linearidades 
quase não distorce. 

Note-se que o sistema XPSK [9] equivale ao 
caso particular em que A,=a; se Aç=]«=1, temos 
o sistema IJF-OQPSK [7], que aliás é também enqua- 
drável no formato OQAM. 
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5. Propriedades úteis dos sinais EDB-OQPSK 


É possível desenvolver sistemas baseados no for- 
mato EDB-OQPSK que realizem um compromisso 
vantajoso entre eficiência de potência e eficiência 
espectral, nas condições correntes da CDVR. Entre 
as referidas condições, são de salientar o uso de 
amplificadores não-lineares e a ocorrência de inter- 
ferências entre canais adjacentes, que podem causar 
uma degradação sensível da qualidade da transmis- 
são se o espectro dos sinais transmitidos for dema- 
siado largo face ao espaçamento entre portadoras. 

Uma das propriedades vantajosas da técnica que 
propomos é o facto de ela permitir, na emissão, a 
síntese directa de sinais que, além de possuírem um 
espectro bastante compacto, apresentam uma envol- 
vente constante ou quase constante — o que significa 
que após uma eficiente amplificação de potência, 
envolvendo um amplificador não-linear, o espectro 
não sofre alteração significativa. 

Além do mais, diversos sinais EDB-OQPSK com 
boa eficiência espectral permitem também uma efi- 
ciência de potência próxima da que caracteriza o 
sistema MSK ideal, padrão de comparação usual em 
CDVR: pode conseguir-se uma degradação inferior 
a 1 dB, relativamente ao padrão referido, transmi- 
tindo sinais EDB-OQPSK com ocupação espectral 
bastante mais reduzida [12]. 


gm, + +—— (8) 
psDO E (A) 
(dB) Mo 

-— 10 NO 


— mm SS ib) 


—— SC le | 


é di (dy 


Mp 
ao 


0 "05 1 |-4T 15 


o os 1 [t=4]T 15 


Fig. 4 — Densidades espectrais de potência (PSD) à entrada 
(A) e à saída (B) de um limitador ideal 
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Alguns espectros de potência de sinais EDB- 
“OQPSK são apresentados. Os espectros foram cal- 
culados mediante a adaptação de um método bem 
conhecido [11] e que se revelou bastante eficiente 
quando aplicado a esta classe de sinais [12]. 


Os espectros calculados (ver Fig. 4(A)) dizem 
respeito a: 


(a) MSK ideal 
(b) IJF-OQPSK 


(c) CPM binário com h=1/2 e sendo g(t) o 
impulso 3RC referido em [2] 


e ainda a sinais EDB-OQPSK caracterizados por (6), 
com 


(d) Aç=e=1/VT 
(d;)) A,=0.8,4=0.6 
(d) A,=0.85, 4=0.53 


Usámos o mesmo método para os sinais à saida 
de um limitador ideal (gm(.)=0, fw(.)= Constante), 
quando os seis sinais atrás indicados são aplicados na 
sua entrada (ver Fig. 4-(B)). 

O facto de os sinais na saída serem também EDB- 
“-OQPSK — por gwni(.) ser constante — facilitou natu- 
ralmente a tarefa. A escolha de um limitador ideal 
deve-se à simplificidade das funções que o caracteri- 
zam e também ao facto de a sua acção sobre os sinais 
se assemelhar à dos usuais amplificadores de potên- 
cia quando operam na saturação. 

Notar que só nos casos (a) e (c) os espectros antes 
e depois do limitador coincidem rigorosamente: isto 
deve-se ao facto de os sinais de entrada exibirem a 
envolvente constante que é própria do formato CPM 
c de não sofrerem portanto qualquer distorção. No 
caso (b), o atravessamento da não-linearidade de- 
grada moderadamente as características espectrais, O 
que tem a ver com a flutuação apreciável da envol- 
vente do sinal à entrada do limitador (3 dB de flu- 
tuação). Nos casos (d;), (d;)) e (d;) a flutuação da 
envolvente é muito ligeira (0.09 dB, 0.18 dB e 0.25 dB, 
respectivamente), e assim não admira que os espec- 
tros de saída e de entrada sejam muito semelhantes. 

Prova-se que é possível conseguir uma degrada- 
ção da eficiência de potência inferior a 1 dB para os 
casos (d,) e (d;) e ligeiramente superior a 1 dB para 
o caso (d,), relativamente ao «padrão MSK» [12]. 


6. Observações finais 


Foi apresentada neste artigo uma técnica de mo- 
dulação para CDVR, baseada num formato aqui 


designado por EDB-OQPSK. 

Esta técnica revela boas características de flexi- 
bilidade e eficiência, aliadas a custos de realização 
relativamente baixos. Utilizando um codificador digi- 
tal é fácil realizar um emissor programável que per- 
mitiria seleccionar qualquer conjunto de funções 
geradoras e qualquer ritmo de transmissão dentro de 
uma vasta gama. Há ainda a possibilidade de usar 
receptores bastante simples, realizados com desmo- 
duladores convencionais do tipo OQPSK: a detecção 
faz-se símbolo a símbolo, por filtragem/amostragem. 
Estes receptores podem ter, para cada código EDB- 
“-OQPSK, um desempenho próximo do desempe- 
nho óptimo que o código permite, bastando para 
tanto que os filtros de detecção sejam adequada- 
mente escolhidos. Um receptor baseado no algoritmo 
de Viterbi pode aproximar ainda mais o desempenho 
óptimo, mas à custa de uma complexidade muito 
superior. 

O formato proposto não só sugere, como vimos, 
uma realização flexível e eficiente do sistema, como 
também possibilita uma análise simples de sistemas 
operando com sinais desse formato: o cálculo dos 
espectros de potência apontado no texto é apenas 
um exemplo. 

A terminar, convém referir que o método de 
codificação /modulação abordado neste artigo pode ser 
generalizado, apresentando no caso geral 4N funções 
geradoras [12]. 
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THE BREWSTER'S ANGLE 
FOR PERPENDICULAR POLARIZA IION 


M. DE ABREU FARO * 
A. S. C. FERNANDES * 


ABSTRACT 


The conditions for total transmission between two 
different mediums (Brewster's Angle) are presented 
for perpendicular polarization of a uniform plane 
electromagnetic wave. An example involving lossy 
mediums also shows a pseudo Brewster's Angle. 


RESUMO 


Apresentam-se as condições de transmissão total 
entre dois meios (Ângulo de Brewster) para uma 
onda electromagnética, plana e uniforme e polari- 
zada horizontalmente. O pseudo Ângulo de Brewster 
é também referido num exemplo com meios com 
perdas. 


The British physicist Sir David Brewster (1781- 
-1868) dedicated most of his life to the study of 
optics, publishing more than 300 papers in his last 
60 years. The Treatise on Optics, published in 1831, 
mentioned his discovery, in 1815, that there was an 
angle of incidence for which the light was totally 
transmitted to the medium that was illuminated 
(1) (2). 

This phenomenon was incorrectly explained for 
a long time on the basis that the electrons excited in 
medium 2 do not radiate in their direction of oscilla- 


* Do Instituto Superior Técnico da U.T.L. e do Centro de 
Electrodinâmica do IN... 
Original recebido para publicação em 15-1-81. 
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tion, which is the direction of the missing reflected 
ray. This argument was destroyed when observing the 
same phenomenon with vacuum as the illuminated 
medium. It is not compatible, either, with the exis- 
tance of a Brewster's Angle for perpendicular polar- 
ized incident waves. 

The explanation of this phenomenon can be sim- 
ply stated involving the matching of the surface 
impedances of both sides of the separation surface. 
Althougth the two different media have different 
characteristic impedances, the surface impedances de- 
pend on the incident and refracted angles, such that 
there is one for which they match, which is the 
Brewster's Angle. 

Considering the theory of the propagation of 
electromagnetic waves, the formal explanation of the 
angle of total transmission (Brewster's Angle) is 
widely presented in the literature. However, only the 
situation of a parallel polarized incident wave is 
normally mentioned. We will briefly show that there 
is also a Brewster's Angle for a perpendicular polar- 
ized incident wave, which is a subject that the authors 
could not find correctly stated in the available inter- 
national literature. 


As shown on Fig. 1, two semi-infinite media are 
separated by the xy plane of a rectangular coor- 
dinate system. The incident uniform plane wave 
has the electric field vector linearly polarized nor- 
mally to the plane of incidence yz, such that it has 
only the x component, E;,. The angle of incidence 
is 6, The media 1 and 2 are described by their elec- 
tromagnetic parameters o-conductivity, e-permittivity 
and u-permeability such that &, = €2/e; is the relative 


El 
e 
Medium 2 a 


[a 


Fig. 1 — Reflection and refraction of a perpendicularly polar- 
ized plane wave. 


permittivity and ty=t,/t, is the relative permeabil- 
ity. The index of refraction of medium 2 relative to 
medium 1 is 


Dz juto (o; + jwe) 
h=>— = (/—————— (1) 
E jom(o, + we) 


and w is the angular frequency. 


The maximum transmission or minimum reflec- 
tion can be calculated from the reflection factor 
Ko for perpendicular polarization, defined as follows 
on the surface of separation: 


E. Ho cos 8, — ni nZ — gin?ô, 
Ro. = a (2) 
Ha cos 8: + y/n2, — sin?8, 


For two dielectric media, expression (2) has a 
zero for the following incident angle 


Fig. 2— €,,, ly space. The shadowed area corresponds to 
the condition of existence of a Brewster's Angle for perpen- 
dicular polarization, for ty, < 1. 


g (V AT ata ) (3) 
= cos! er ra 
| l—u, 


which is the Brewster's Angle for perpendicular polar- 
ization. Analysing this function, it was concluded 
that there is a Brewster's Angle for perpendicular 
polarization when the following conditions exist 
simultaneously for the media present 


Ho << E 
<1/e (4) 


in the basis of ; <1. 

Figure 2 shows the corresponding area of the &,, 
>, space, for the case of  < 1. 

For the hypothesis of ta > 1, as obvious, we 
should have 2, > €&, and uy > 1/6». 
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Fig. 3 Modulus and phase of the reflection factor for 
perpendicular polarization. =, =28, 4=01, 0,=0, 
o, = 105 m'. 
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When the media are not purely dielectric there 
is, as for the parallel polarization, an angle of mini- 
mum reflection called the pseudo Brewster's Angle. 

For the media characterized by the point A on 
Fig. 2, &,=2.8 and w,=0.1, Fig. 3 shows the modulus 
and phase of the reflection factor for two frequen- 
cies. This example considered c,=0 and ,=107S m”. 
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AN OPTIMAL POLICY FOR CONTROLLING 
NON-INTERACTIVE FINITE SEIS 
OF ELEMENTARY PROCESSES 


JOSÉ QUINTINO ROGADO, 


Full Professor of Instituto Superior Técnico. CVRM (INIC) 


ABSTRACT 


An optimization method to deal with finite sets of nonin- 
teractive elementary processes is presented. 

The objective function is the rdtio of the intensive varia- 
ble to one extensive variable. 

The method can be used with any physical discrete 
systems, where the extensive state variables are non increa- 
sing functions of the intensive state variable. 

The optimazing path is directely difined avoiding the 
heavy algorithms usualy used in non-linear programing. It 
is recomended for situations where very numerous data 
should be processed, as it riquires reduced computer me- 
mory capacity. 


SUMÁRIO 


Descreve-se método para optimizar conjunto finito de 
processos elementares não interactivos onde a função ob- 
jectivo que se oferece é o ratio entre a varidvel de estado 
intensiva e uma varidvel de estado extensiva. 

O método, aplicável a diversos sistemas físicos discretos 
ou discretizáveis onde as variáveis de estado extensivas se- 
jam funções não crescentes da varidvel de estado intensiva, 
recorre à definição directa da trajectória optimizante dis- 
pensando os pesados algoritmos habitualmente usados em 
programação não linear. É recomendado para situações 
onde o número de dados é muito abundante, pois exige 
reduzida capacidade de memória computacional, 

|. Let us consider a discrete system S formed by a finite 
set of n non-interactive sub-systems, fed by n elementary 
sources, each one allocated to a sub-system. 

At any functioning point, S continuously yields an out- 
put if it is continuously fed by an input source made up by 
the weighed union of the elementary sources. The process 
carried out on S will be called S process; the ones developed 
in sub-systems s;, elementary processes s;. 

In order to be concerned to actual physical situations, it 
will be assumed that the imput source will be finite. Then, 
the output product will also be finite. 

In many instances the control of the system output can 
be made using three positive non-independent macroscopic 


state variables. Two extensive affine ones, X and Y, and an 
intensive one, Z among which there exists the fundamental 
relationship. 
X= ZY A TZC R. (1) 
Yet, at the elementary level the process carried out by s; 
sub-system is controlable by three positive non-indepen- 
dent state variables. The extensive ones, x; and y;, and the 
intensive one, z; (:x;, Y;, 40 R$). 
In paralel with (1), there are n equivalent relationships, 
once the “qualities” of each elementary source are fixed: 
X=2ZY% d1eftl,2,..,n) (2) 
Note that elementary sources can belong to a time sequen- 
tial stream, i.e., a time variable source forming a time se- 
ries. In this case, the elementary processes actuate sequen- 
tially over the source time series. 


In any case, as the elementary processes are not interac- 
ting, the following stays: 
X=hix= Aly; z 
= Aq je [1,2,..n) 
(3) 
dA vi my 


n ê 
Ze za s ai=1 
A! y; | = 


Where, A! are weighting factors to deal with the allocation 
of the elementary sources into the overall source. 

In the sequel, it will be assumed that whatever 115, both x; 
and y; are functions of z;. This is, s; processes are such that 
z; can be taken as the independent state variables. 

Taking into account (3), it is obvious that the S process 
can be defined using as state variable the vector 
z (z1, Ed em Zn) 1 £e RE : 


X = X (Z) 
Y=r(Z) (4) 
Z=Z(D 
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2. Let assume that the state of any s; process is defined 
overa finitesetof points zy eR : ke [1,...,n;j and that 
(zw) el, being || C R 4 bound intervals. 

Then, the evolution of S can only be known over discrete 
points of space state corresponding to U [z;). 

For simplicity sake, the initial state of every s; is defined 
by Z; = Zijo, Such that: 


a) Zjy is an attainable state 


b) Y ke [l....mi 
any attainable state, 


Zio É Zik Where z; is 


Very often, the physical meaning of state variables leads 
to the situation where both x; and y; are non-increasing 
functions of z,. 

In this case, the initial state variables x;, and y;o are the 
upper limits of the respective state variables. 

Morever, the initial state of S is defined by Z, 

such that: 


a) Z, is an attainable state 
b) W attainable Z:Z, 5 Z 


Considering (3) and assuming that x, and y, are non-in- 
creasing functions of z;, the initial state corresponds to 
X = sup [X) and Y = sup [Y). 

If S isina certain state Z and is transfered to another 
state, let assume that two kinds of transitions (trajectories) 
in state space are possible. The ascending transitions, when 
Z increases (and X and Y do not increase); the descending 
ones, when the contrary occurs. 

Itis obvious that, if S is in the initial state, it can only be 
changed along an ascending transition. 


3. Among all the possible transitions of S, starting from a 
fixed state, the one ensuring the following maximum local 
condition 


NL 
max AX (5) 


is called an optimal transition. 


The physical meaning of an optimal transition is the fol- 
lowing;: the optimal transition ensures the maximum increa- 
sing (decreasing) in the intensive state variable Z with the 
minimum decreasing (increasing) of extensive variable X. 

Of course, the optimal transition can be either ascending 
or descendineg. 

In order to build up optimal transitions, lets consider the 
elementary transition of S from states Z, to Z,, only caused 
by the transition of one elementary process, say s,, from s,| 
LO Spo. 

The extensive state variables of s. will, thus, undergo the 
following variations: 


012 Xe = Xez — Xel 


012 Ve = Yo — Yel 


By (3) the corresponding variations of S extensive state 
variables are: 


A 12 X = A" (x — Xe) = Aº 012 Xe 
Ao Y A (vez — Ye) = Aº 612 Ye 


II 


Theorem I; An elementary transition of s; from state Z;x 
to state Z;. is an ascending transition of the process S from 
states Z, to Z. if and only if: 


(zx — Zy) Yi < Ci — Ly) Vie (6) 


Proof: According with (1) and (3), the elementary transi- 
tion does not cause any transition on S if 


Ar Yr t Aj (Vie — Yir) E 


[q 
ae 
tas 


Ar YrZr + Ai (Vie Zie — Yik Zik) 


| 
DE neto 


which means that 


Z, Ai (Vie — yik) = Ai (ie Zie — Yik Zik) 
or 


Vie (2x — Zk) = Yie (Ze — Z,) 


Then, if (6) helds, there is a transition in S from Z, to Z, 
andZ, > Li. 

Of course, a similar theorem can be demonstrated for the 
descending transition “mutatis mutandi”, 

Lemma: If x, = O, then y, = O. This lemma results 
imediately from (2) and from the fact that z; is bound. 

Theorem II: If x, are non-increasing functions of z; and 
If x; (y;) pass in the origin, the functions y; (x;) are convex 
functions. 

Proof: If x; (z;) is non-increasing, then. 


z; (Ax;)) = 2; (x;) 


or 
y; (À X;) = MH (x;) dr é Ú 
À x; X; 
or 
vi lhx;) ES Ay; (x) 


which is the condition of convexity for a function passing in 
the origin. 


Corollary: If y; (x;) passes in the origin and if x; (z;) is 
a non-increasing function, y; (2;) is a non-increasing func- 
tion. 

Proof: By Theorem Il, y; (x;) are convexes, morever they 
pass in the origin. Then, they are non-decreasing; x; (z;)) are 
non-increasing. Thus. 


Yi (Xi; (C)) = yy (2) 
are non-increasing functions. 


In the sequel, it will be assumed that y; (x;) pass the origin, 
in order to be in a actual physical situation. Yet, x; (z;) and 
y; (z;) are non-increasing functions, and y; (x;) are con- 
vex ones. 


If, in any elementary process s;, an ordered relationship 
1s built up. 


EI E tp ES ve ME Tate [Za] 


and if s;is in the state z,,, obviously, the transition to Z; + | is 
an ascending one. 


Let consider the above mentioned one-step transition 
operator E= | + À and let define the ratio R as follows: 


R (85) = AttttMi MG Oh rt Yin) Yirtt o Vir 


O rctiX AMX rtt Xi) Mir+ — Kir 


Of course, E causes ascending transitions on process s; 
É (6) is fulfilled. The result of the E operator can be 
mesured by the ratio of increasing of Z in terms of decrea- 
sing of X. 

ltis worthnote that the E operator actuate only ordered 
sets [z;), as stated above. 


4. Lets assume that the process S is in its initial state z, and 
that all one-step elementary transitions carried out by E 
operator are considered. 

Theorem III: Ifthe finite set of one-step transitions satis- 
fying (6) is arranged in a decreasing ordered sequente of 
ratio R, the first possible elementary maximal increasing 
transition of process S is the one corresponding to the first 
term of the above mentioned decreasing ordered sequence. 
Yetthe X (Y)thus built up is the upper convex ofall X (Y) 
functions. 


Proof: This optimal ascending policy results immedia- 
tely from Theorems I, and II and from its corollar and from 


the obvious relationship between. SE and A 
that Z and Y belong to the L classe functions 


assuming 


Of course, the optimality principle of Bellman is applied 
to this policy principle. This means that the policy is valid 
for any state of the S process giving rise to a recurrent 
policy. 

An optimal descending policy can be demonstrated in 
similar terms “mutatis mutandi”, 

5. The presented optimal policy principle allows to build 
up he optimal trajectory of S process avoiding either the 
use of strong combinatory algorithms or the combersome 
non-linear optimization methods, when the envolved func- 
tions are not affine. 

As a matter of fact, it delivers an easy way for construc- 
ting the upper convex of all X (Y) functions with relatively 
small computer memory consumption. 
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CONDUTIVIDADES IÓNICA E ELECTRÔNICA 


DUM ÓXIDO PURO 


MOy COM DESORDEM DE SCHOTTKY * 


C. A. €. SEQUEIRA 


Laboratório de Electroquímica, Instituto Superior Técnico, Lisboa. 


RESUMO 


Usando o método de aproximação de Brouwer, 
calculam-se as concentrações de defeitos para o caso 
dum óxido puro MO, com desordem de Schóttky em 
função da pressão parcial do oxigénio na fase gasosa 
em equilíbrio com o óxido. As expressões matemá- 
ticas das condutividades iónica e electrónica, bem 
como do número de transporte, que se podem dedu- 
zir das concentrações são também apresentadas. 


ABSTRACT 


Appropriate defect reactions are formulated and 
equilibrium defect concentrations for a pure metal 
oxide MO, with Schônky disorder in equilibrium 
with a surrounding oxygen gas phase are calculated 
using the Brouwer's method of approximation. Mathe- 
matical expressions for the conductivities and trans- 
port numbers of the predominating ionic and elec- 
tronic charge carriers are also deduced. 


INTRODUÇÃO 


Seguindo o método de aproximação de Brouwer, 
o qual constitui uma ferramenta típica na análise 
prática de problemas de química do estado sólido [1], 
apresenta-se neste artigo o cálculo das concentrações 


de defeitos para o caso dum óxido puro MO, com 
desordem de Schottky em função da pressão parcial 
de oxigénio na fase gasosa em equilíbrio com o 
óxido. Estabelecem-se também as expressões das con- 
dutividades, bem como dos números de transporte, 
que se derivam prontamente das relações das con- 
centrações. Trata-se duma situação já analisada por 
Tallan e colaboradores [2] para compostos de fór- 
mula geral M;Xs, MX, M;X, M;X; e MX>, mas que 
no presente trabalho é generalizada a óxidos do tipo 
MO, com y = 1/2, 1, 3/2, 2, etc. 

Simplifica-se este estudo admitindo que o óxido 
MO, é essencialmente iónico, possuindo os seus defei- 
tos pontuais cargas efectivas normais (isto é, não se 
considera a ionização desses defeitos) e que não existe 
essociação de defeitos. Trata-se duma situação bas- 
tante frequente em estudos de sistemas do tipo óxido- 
-gás a altas temperaturas. 


Para facilitar a comparação com outros trabalhos 
publicados sobre assuntos análogos, usa-se a notação 
atómica de Krôger e Vink [3]. 


EXPRESSÕES DAS CONCENTRAÇÕES 


O sistema MO, — O, a alta temperatura, no caso 
em consideração em que, portanto, a desordem pre- 
dominante é uma desordem de Schóttky, é sede dos 
seguintes. equilíbrios: 


* Comunicação ao IV Simpósio Brasileiro de Electroquímica e Electroanalítica, São Carlos, Brasil, Abril de 1984. 


original recebido para publicação em 19/6/84 
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O=e+h' (1) 
y 

3 0,=Vi +2yh' +yo, (2) 
O=V3 +yv: (3) 


em que O, representa um átomo de O num lugar 
normal da rede cristalina do óxido MO,, VY” repre- 
senta uma lacuna de M dotada de carga efectiva —2y, 
Vo representa uma lacuna de oxigénio dotada da 


carga efectiva +2, e' e h' representam electrões livres 
e buracos de electrões com as cargas efectivas —1 e 
+1, respectivamente, e O designa o cristal perfeito 
MO,. Designando por K,, K;, K; as constantes de 
equilíbrio correspondentes respectivamente às rela- 
ções (1), (2) e (3), por P(O,) a pressão de oxigénio e 
por ne p as concentrações de electrões livres e bu- 
racos de electrões, podem deduzir-se as três relações 
seguintes: 


n.p=kK, (4) 
| Vê |. p= MO) (5) 
| de | | Va | Y=K, (6) 


em que as concentrações dos diversos elementos 
estruturais se exprimem em fracções de lugares da 
rede cristalina. As relações básicas acabadas de escre- 
ver, pode juntar-se uma 4.º relação tradutora da 
electroneutralidade do sistema, a qual se escreve: 


2 |Vy|+n=2 | Vo | +p (7) 


Resulta então um sistema de quatro equações a qua- 
tro incógnitas — as concentrações n,p, |V”"'| e |[Vo| — 
que são calculáveis a uma dada temperatura T, em 
função de P(O;) e das constantes K,, K; e Ks. De 
acordo com o método de aproximação de Brouwer 
[4] o cálculo pode simplificar-se desde que se consi- 
derem domínios de P(O;) tais que em cada membro 
da relação (7) só uma das concentrações seja 
predominante. Obtêm-se assim três simplificações 
possíveis: 


l. Domínio de P(O,) elevadas — neste caso o 
equilíbrio (2) está muito deslocado para a 
direita, resultando p>n e [VR |» | Voe A 
equação de neutralidade reduz-se a: 


2y | va | =P 
2. Domínio de P(O;) fracas — agora tem-se p <n 


e [VW |«|V9' de modo que a relação de 
electroneutralidade se reduz a: 


2 |Vo |=n 


3. Domínio de P(O;) intermédias — neste caso é 
preciso distinguir duas situações: 


(a) Se K, > K,;, n > Ed ep > dim é de 
modo que a equação de neutralidade se reduz 
a: 


n=>p 


(b) Se K, « K,n «< | V3"] ep < |Vi| vindo 
para a equação de neutralidade a relação: 


PIYal= vg 


No domínio de P(O,) elevadas predominam as lacunas 
de M e os buracos de electrões complementares, de 
modo que os óxidos se comportam como não este- 
quiométricos com déficit de metal (M,« O). Como 
exemplos de óxidos que apresentam déficit de metal 
e apreciável condutividade iónica a pressões parciais 
de oxigénio elevadas podem referir-se o V;O;, Cu,0, 
NiO, etc. Note-se que o NiO é a roximadamente este- 
quiométrico a fracas pressões parciais de oxigénio o 
que está de acordo com o facto de que num mesmo 
óxido se podem observar diferentes «situações defei- 
tuosas», consoante a pressão de oxigénio da atmos- 
fera envolvente do óxido. 

No domínio de P(O;) fracas predominam as la- 
cunas de oxigénio e os electrões complementares, 
havendo um afastamento da estequiometria que se 
traduz pelo aparecimento de redes M,.x O, com 
déficit de oxigénio. No caso destes compostos (con- 
dutores mistos), observa-se uma transição da sua 


condutividade electrónica do tipo p a elevadas P(O,) 
para uma condutividade do tipo n a baixas P(O;). 
Esta situação é típica do Ta,O,, V,Os, ZrO;, etc. 

No domínio de P(O;) intermédias, para o qual se 


tem K, > K;, | V2”| aumenta e | V:*| diminui à me- 
dida que P(O;) aumenta, enquanto que n e p perma- 
necem constantes. Trata-se duma situação particular- 
mente importante nos sistemas de óxidos em que a 
ionização intrínseca é elevada sendo mais importante 
que os desvios à estequiometria, predominando as 
concentrações dos portadores electrónicos na rede 
cristalina (ex. CaO — ZrO, a altas temperaturas). 


Finalmente, no domínio para o qual se tem 
K, < K;, n decresce e p aumenta à medida que P(O,) 
aumenta, enquanto que |V) | e | VX”| permanecem 
constantes. E uma situação típica de electrólitos sóli- 
dos em que predominam os defeitos de Schôttky. 
Como exemplo, pode referir-se a tória pura que exibe 
condutividade iónica predominante a pressões parciais 
de oxigénio muito reduzidas. 
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No quadro I mostram-se os resultados obtidos por 
resolução do sistema de equações resultantes para os 
diversos domínios de pressão parcial do oxigénio em 
que se desdobrou a condição de electroneutralidade 
para o caso em estudo. 


EXPRESSÕES DAS CONDUTIVIDADES 


A condutividade eléctrica dum cristal deve-se à 
presença de defeitos pontuais carregados na rede 
cristalina, defeitos esses que podem ser de natureza 


electrónica (e” ou h') ou atómica (p. ex, VE E Vad 
A condutividade eléctrica devida a uma partícula 
carregada i está ligada à mobilidade desta partícula 
pela relação: c = c; z, e v;, onde c, é a concentração, 
z; a valência, v;, a mobilidade da edema jecéa 
carga do electrão. Uma vez que para os óxidos MO, 
que estamos a tratar, a corrente iónica é transportada 
por lacunas de O e de M, sendo os electrões livres e 
os buracos de electrões os responsáveis pela corrente 
electrónica, pode escrever-se: 


Go =(VE)+HoVI) o =o(e)+o(h) 


em que (VI, S(V.), o(e!) e o(h') são as conduti- 
vidades das partículas V$”, Vyp e e h', respectiva- 


mente. 


Representando as mobilidades dos electrões, dos 
buracos de electrões, das lacunas catiónicas e anió- 


nicas por v(e”), v(h'), v(V59 e vCV es respectivamente, 
e atendendo à forma como estamos a definir as 
concentrações de defeitos, podemos então exprimir 
as condutividades iónica e electrónica, cj € c«, pelas 
relações que constam da Tabela I, em que K,, Ks, 
K, e K, são constantes dadas pelos valores tabelados, 
para as diferentes aproximações da condição de neu- 
tralidade Se a mobilidade dum defeito é indepen- 
dente da sua concentração, a condutividade resul- 
tante é directamente proporcional a essa concentração. 
Esta proporcionalidade verifica-se geralmente no caso 
dos óxidos puros desde que a concentração em 
defeitos carregados seja reduzida. Então, as equações 
da Tabela I permitem exprimir as leis de variação 
das condutividades iónica e electrónica em função 
da pressão parcial de oxigénio em equilíbrio com o 
óxido, tendo as suas variações logarítmicas um anda- 
mento aproximadamente rectilíneo e análogo ao que 
se obtém nas mesmas condições para as concentrações 
de defeitos em função de P(O,). 


EXPRESSÕES DOS NUMEROS DE TRANS- 
PORTE 


Num composto inorgânico de condução mista, os 
números de transporte das espécies iónicas e electró- 


Tabela | Expressões deduzidas para as várias aproximações da condição de neutralidade. 


Bis at? Xe H2y 213 P( 0,) 


” Kika 


2 "[K, 


-113 near 1223, 


(K, RS 2 213 
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116 -12/y+], 
y is 


N2y M2yly Ala 7 
2 K ho) | 


244 2 viva 2(yk, sá Prlvj) 
a vet) 


nf? vth') 


Ez 
tK 


nicas, ti; € ta, estão ligados às condutividades pelas 
relações: 
tio = cjo/€ ta = 0/0 

uma vez que c=csotocy=0(toTt ta) com to + 
+ ta = 1. Por conseguinte, as expressões matemáticas 
do número de transporte podem deduzir-se facilmente 
a partir das relações das condutividades electrónica 
e iónica, acabadas de estabelecer. Tem-se, para ns 
números de transporte: 


l. Caso em que 2y |VY'| =P 
E K, P(0,;)” + Ks o 
to (Rs + K9 PO)? + Ks + K (O 
oe Ks P(O;)'* + K; P(0))” 
e TT 


(K, + K;) P(O))* + Ks + Ki P(O)* 


2. Caso em que 2|Vo/=n 


DO Mk POrr 
o RA + (Ks + Ky P(Oj IP + K, P(O)IP 


Ks P(0;)'” + K; P(O;)'” 


“Ra + (Ks + K9 P(O91P + K, PO)? 


3. Caso em que n=p 


K, P(O;)'? + Ks P(O)) 1º 


o TA P(O,VR + Ks (OP + K, + K; 
K +K; 
ta = 


K, P(O,Y2 + Ks P(O)) 12 + Ks + K; 


4. Caso em que y |V3'| = |V5'| 


o K, + Ks 
o = K, + Ks + Ks P(O;))!" + K; P(O)!! 
o K, P(0;) 4 + K; P(O;)!* 
let = K, + Ks + K, P(O)1! + K,; P(O) 


As expressões aqui deduzidas, embora constituam 
uma solução aproximada, permitem determinar, duma 
maneira simples e bastante rigorosa, propriedades de 
óxidos MO, e processos neles ocorrentes como os 
de transferência de massa, reacções óxido-gás, etc., 
facilitando a discussão prática de inúmeros problemas 
de química do estado sólido, como refere Krôger [1]. 
No caso particular que tratámos de que já se fez uma 
comunicação [5], os cálculos apresentados constituem 
um primeiro approach a um estudo mais geral sobre 
as propriedades eléctricas de sólidos a alta tempera- 
tura. 
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IWO STUBS ANDA/4 LINE IMPEDANCE 
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RESUMO 


Nesta nota descreve-se um procedimento para con- 
seguir a adaptação de impedâncias com a montagem 
clássica constituida por dois «stubs» e uma linha 
de h/4 de comprimento. Verifica-se por meio de 
exemplos que a convergência é rápida em quase 
todos os casos e que a simplicidade do processo per- 
mite que seja utilizado por operadores sem qualifi- 
cações especiais. 


ABSTRACT 


In this note a simple procedure to match a load 
to a transmission line with the two stubs and A/4 
line set up is described. This procedure is normally 
quickly convergent and can be carried even by untrai- 
ned operators. 


1. INTRODUCTION 


The use of a set up consisting of two stubs and 
a quarter wavelength fixed line to match a load to 
a transmission line at a given frequency is well known 
[1]. The appropriate lengths of the stubs are easily 
calculated as long as the admittance of the load is 
known and in many instances the implementation of 
the set up presents no difficulties. 

However there are cases in which this is not true. 
One such case arises when coaxial air-filled lines 
are used. In these lines the central conductor is sup- 


Original recebido para publicação em 11/12/84. 
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ported by dielectric beads with the effect that elec- 
trical length is not simply converted to mechanical 
length but depends on the number of beads (of uncer- 
tain permittivity) inserted in the line and in the 
connecting Ts. Thus same amount of fiddling is 
needed to translate one length into the other, resul- 
ting in time consuming and in the possibility of 
making mistakes when determining the length of the 
stubs and of the quarter wavelength line itself. 

Most often, as a result, when the set up is tested, 
the matching is very far from being reached and the 
final position must still be obtained by trial and error. 


We have had this experience repeatedly in the 
laboratory classes where often the students were una- 
ble to obtain the match, with the ensuing frustration. 
Thus we were led to devise a simple “non think” 
matching procedure entailing only SWR measure- 
ments and without requiring any previous settings 
of stub length. Moreover in most cases it converges 
quickly to a predetermined value of the SWR. 

In spite of its extreme simplicity this procedure 
has not been mentioned and much less its conver- 
gence studied before, as far as the authors are aware. 
As it may be useful in situations not only as the 
one mentioned above but also where non technical 
staff is used, it was decided to publish this short note. 


2. DESCRIPTION OF THE PROCEDURE 


The procedure consists simply in adjusting, in 
successivc steps, each stub in alternance in such a 
way as to obtain a minimum SWR each time. 

As a preliminary step and just for the sake of 


the argument, we assume that, stub A (fig. 1) is 
adjusted so as to reduce to zero its input admittance 
(open-circuit). Then the first step consists of varying 
stub B until minimum SWR is measured in the 
transmission line. In the second step it is the length 
of stub A that is changed until a new minimum 
SWR ensues. In the third step it is again the turn 
of stub B to be changed, and so on. 


Tranambanios lina 


tdi Led (b) (A) Load 


In figures 2 several examples of matching by this 
procedure are shown. In figure 2a the paths P, repre- 
sent the locus of the admittance at terminal plane 
(d) (fig. 1) when the length of stub B is changed 
and the paths Q, represent the locus of the admit- 
tance at the same terminal plane when the length of 
stub A is changed. Path R represents the locus of 
the admittance seen at terminal plane (b) when the 
length of sutb A is changed and path S is its image 
at terminal plane (c). In the remaining figures these 
notations have been omitted. 

It can be seen that this procedure is convergent 
to the center of the Smith chart (i.e, matching): by 
varying stub A each time in such a way as to mini- 
mize the SWR on plane (d) the point representing 
the admittance at plane (d) falls on a constant con- 
ductance curve which every time is nearer to the curve 
of conductance g = 1 (except, of course, when this 
point is on the curve g = 1 itself). 

As can be seen the procedure converges quickly 
in the sense that few steps are needed to transform 
an original load with a large SWR to a final load 


Fig. 2a 


with small SWR. On the other hand the procedure 
is experimentally straighforward: the minimum SWR 
to be obtained in each step is easily identified 
and even if the exact location of the minimum 
SWR may be difficult to obtain (as is the case, 
for example, in step number 2 of figure 2c) the 
convergence is not greatly affected. 

Now it should be obvious that the initial value 
of stub A is of no consequence as its initial admit- 
tance can always be assumed to be included in the 
value of the load admittance. 


3. SPECIAL CASES 


Some special cases need attention, One of them 
arises when the load admittance is real. In the gene- 
ral case, as can be noticed from the previous exam- 
ples, step 1 brings the admittance at terminal plane 
(d) to a real value. Thus in the present case this 
step is omitted. But then, as can be easily verified, 
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the change of stub A in any of the two possible 
directions has the effect of increasing the SWR lea- 
ving the idea that matching is impossible. There is 
neverthless a simple remedy that consists in giving 
to the admittance of the load a reactive component 
using stub A for this purpose. Then step 1 follows 
as for any other load. 

When the load is almost real convergence is still 
possible but becames very slow: see fig. 3. The 
remedy is again to give a larger reactive part to the 
load with the help of stub A. 


FIG: 


In the conditions described in the introduction it 
happens quite often that the fixed line is electrically 
different from A/4. As is known, matching with the 
two stubs and fixed length line set up is still possible, 
but as can be seen from the comparison of figures 4 
drawn for a À/8 fixed line with figures 2 drawn for 
a A/4 fixed line, the convergence becomes slower. 
Thus as far as convergence is concerned, it is advan- 
tageous to use a, as near as possible, quarter wave- 
length fixed line. 


FÃS 
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Fig. 


Fig 


4. THE SECOND SOLUTION 


It is well known that, for each load admittance 
that can be matched, this set up has two solutions in 
the sense that matching is obtained for two different 
pairs of values of the admittance of the stubs [1]. 
However the procedure described before converges 
without ambiguity to only one of these two solutions. 

To find the other solution a divergent step must 
be introduced somewhere. An example is shown in 
figure 5 in which the admittance to be matched is 
the same as in figure 2a. Now, it is the admittance 
of stub B that is assumed to be initially zero and a 
first convergent step is taken with stub A. But then 
the second step (with stub B) is divergent and stops 
at an arbitrary value of stub B admittance. Onwards 
the other steps always converge to the central point 
of the Smith chart. 

Clearly this procedure looses the deterministic 
characteristic of the previous one: in step 2 there 


Fig. 5 


is no way of choosing the sense in which stub B 
must be changed to lead to a converging situation 
(either to points 2 or 2D in the figure) and there is 
no rule saying when to stop changing the length of 
this stub. Thus it does not seem to be of much 
practical interest. 


5. CONCLUSIONS 


A simple matching procedure for a two stub and 
quarter wavelength fixed line has been described. It 
is seen, through various examples, that in most cases 
the procedure converges quickly when the fixed 
length line is equal or near to À/4 and more slowly 
when it departs from this value. The procedure has 
the following advantages: (i) it does not need any 
previous calculations or initial settings (except for 
one of the stubs which must be at open circuit or 
near): as a consequence it does not entail the trans- 
lation of electrical lengths into mechanical ones and 
the exact length of the fixed length line is of secon- 
dary importance; (ii) it is very simple to follow as 
it implies only the measurement of standing wave 
ratios — as a matter of fact it can be executed by 
an operator without special training. 

Several years of experience in the laboratory clas- 
ses have shown that the procedure is effective. 
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THERMAL RECTIFICATION IN EXTENDED 


CAVITIES 


J. M. M. BARATA, 
D. F. G. DURÃO 


Instituto Superior Técnico, Mechanical Engineering Department 


RESUMO 


Um programa numérico que resolve as equações 
da energia, função de corrente e vorticidade foi 
usado para quantificar o efeito do número de Raleigh 
e geometria na rectificação térmica em cavidades. 


ABSTRACT 


À numerical program solving the energy, stream- 
-function and vorticity equations was used to quan- 
tify the effect of Raleigh number and geometry in 
the thermal rectification obtained with extended 
cavities. 


INTRODUCTION 


Available information on the velocity field and 
heat transfer in square and rectangular cavities has 
been lately published most related with insulation and 
solar energy absorvers (Barata, Durão and Reis, 1983; 
Koutsoheras and Charters, 1977; Lopes, 1981; Ozoe, 
Sayama and Churchill, 1974; Quere Humphrey and 
Sherman, 1981; Shiralkar and Tien, 1982). 

The present work is motivated by the desire of 
knowing the capabilities of using cavities to improve 
the heat conservation and insultion of buildings. It is 
pretended to quantify the increase of thermal recti- 
fication in extended cavities as these shown in Fig. 1. 


original recebido para publicação em 20/1/85 
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Barata et al (1983) have already shown that the addi- 
tional spaces acting on the fluid movement and con- 
vection increase the thermal rectification. 

A numerical program solving differential equa- 
tions for energy, stream-function and vorticity was 
used to obtain the velocity and temperature fields 
from which the Nusselt number was computed. Gri- 
des of 32X32 points were used and the conversion 
was tested and optimized considering square cavities 
(Barata et al, 1983). 


FIG. 1 


K 


Next section presents the solved equations and 
the conclusions follow the section were the numerical 
and experimental results are introduced and discussed. 


SOLVED EQUATIONS 


The following non-dimensional Navier-Stokes 
energy and continuity equations were solved: 


Pr OX dy 
+ Ra (cos = sing) 
OX Oy 
oT oT | 
——— + ————. 2T 
E 0X É Oy ii 
a em A 
RA 
dy 
E E qe 
OX 


RESULTS AND DISCUSSION 


The next figures represent a set of results which 
report the main pictures of the detailed investigation 
carried on. Fig. 2 shows for H,=H/2 and L,;=L/5 
the dependence of Nusselt number on Raleigh num- 
ter for Ra between 3X 10º and 2X 10. The results 
are presented for values of B (rotation of the coor- 
dinate system) equal to 90 and 270 degrees which 
corresponds to the inversion of gradients, i.e. the 
substitution of cold surface by hot surface and vice- 


“versa. An example of such conditions is the alterna- 
tion of days and nights and the consequently inver- 
sion of heat fluxes verified in the wall of buildings. 


In can be recognized that the thermal rectifica 
tion is present and the efficiency, defined as 
N N 


Ug=0) "NU cg.2m) 


= 
N 


U cg =20) 


increases with Raleigh number. For Ras X 10º, 
verified in practical situations with low surface tem- 
peratura differences (AT), the efficiency rounds 239%. 
The variation of Nu and n with Ra is: 


(8 = 90) Nu = 0.0402 Raº-3% 
(8 = 270) Nu = 0.0685 Raº317 
n=1-— 1.705 Ra0.973 


The efficiency obtained with Ra = 50000 and 
H,/H, = 1.5 is 19.6%; further investigation demons- 
trated that the maximum efficiency corresponds to 
values ot H, = Ho. Fortunately this geometry is the 
one that can be feasible and easily manufactured. 

Fig. 3.1 and 3.2 show for adiabatic top and 
bottom walls, H,=H,, L,=L/5, Ra=2.7X 10 
and Tho/Tooda = 30/10 (Pr = .71) isolines of non- 
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-(limensionalized stream-function and temperature 
corresponding to 8 = 90º and 270º respectively. 


They visualize the flow for cases of favourable 
and non-favourable heat-transfer. For B = 90º the 
resulting Nu is 4.30 in contrast to the value of 2.95 
obtained for 8 = 270º (mn = 31.4%). This large ther- 
mal rectification is due to the existance, for 8 = 270º, 
of two vortexes which by creating two quasi-inde- 
pendent flows diminuish the mass and heat transfer 
between the two isothermal walls. For the geometry 
of Fig. 4, Table 1 introduces the effect of AT. 


Fig. 4 and Table 2 present the influence of H/L 
indicating the existance of a maximum value of n 
for H/L around 1.5. It should be emphasized that 
the elongation of the cavities also causes an impor- 
tant diminution of the heat transfer even for B=90º. 
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H/L 


TABLE 1 — Effect of àT on m 


Nu (8=90)  Nu(8=270º) 


2.873 2.448 
3.660 2.620 
3.930 2.142 
4.300 2.950 


TABLE 2 — Effect of H/L on m 


Nu (8=90) — Nu(B=270º) 


2.873 2.448 
1.791] 1.401 
965 899 


n 


14.8% 
28.4% 
30.2% 
31.4% 


Table 3 gives the influence of L;,/L in the values 
of efficiency and thermal rectification indicating their 
increase with L;; the heat transfer diminuishes with 
the increase of L,. 


TABLE 3 — Effect of L,/L on n 


L,/L Nu (B=90º) 


Nu (8=270º) n 
O) 3113 Sds O) 
1/5 2.813 2.448 14.8% 
1/4 2.694 2.219 17.6% 
CONCLUSIONS. 


The value of the thermal rectification increases 
with Raleigh number and dimension L,. The varia- 
tion of the thermal rectification with H/L has a 
maximum, which for L, = L/5 is located at H/L 
around 1.5. 
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SOBRE O CONSUMO ENERGÉTICO 
DOS VEÍCULOS AUTOMÓVEIS 


(UM COMENTÁRIO) 


F. ALMEIDA E CASTRO, IST. 


Segundo números divulgados internacionalmente 
[1], Portugal situa-se na área da OCDE como um 
dos países mais consumidores de energia nos trans- 
portes, relativamente ao seu Produto Interno. Por 
singular que pareça, estamos neste domínio muito 
perto do grupo dos mais perdulários (E.U.A., Canadá, 
Irlanda, Austrália) e a considerável distância do 
pelotão dos mais eficientes (Escandinávia, Alemanha 
Federal, Japão); na Europa, além da Irlanda, só a 
Turquia nos tem ultrapassado nesta corrida para a 
penúria. Não é uma posição lisongeira, tratando-se 
do sector mais vulnerável da nossa economia — de- 
pende do petróleo em 99% — e contribuindo em 
forte medida para o nosso endividamento externo. 

Ocorrendo que a maior fatia do consumo de car- 
burantes no transporte cabe ao automóvel — em 
Portugal responde por 85% da energia dispendida 
na movimentação de pessoas — surpreende que às 
questões relacionadas com o rendimento energético 
dos automóveis se conceda, entre nós, bem pouca 
atenção; tirando as ocasionais referências que sur- 
gem na imprensa dita especializada — na sua maioria 
simples transcrições mais ou menos laudatorias de 
textos publicitários de origem estrangeira, sem a dese- 
jada correlação com a especificidade do caso por- 
tuguês — o vazio é regra no que respeita ao conhe- 
cimento, difusão e debate sobre as componentes efec- 
tivas do problema, em especial por parte de quem 
(por formação de base) melhor poderia fazê-lo com 
objectividade. 


Original recebido para publicação em 21/12/84, 
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Vêm estas considerações a propósito da recente 
publicação numa conhecida revista inglesa de divul- 
gação científica dum texto em que a questão do 
consumo de combustível nos automóveis é abordada 
na ótica do utente [2]. Suportando, com algum exa- 
gero, a tese de que o factor aerodinâmico não tem 
a mínima relevância na prática corrente, avança Oo 
autor do artigo uma espécie de directiva global para 
uso dos potenciais interessados; esta contudo, por 
excessiva simplificação e deficiente enunciado, 
arrisca-se a confundir mais que a esclarecer. 


Parecendo de interesse, pelo que já foi dito, desen- 
volver este mote com uma preocupação de maior 
clareza e sem perder de vista as reais perspectivas 
do País em que vivemos, seguidamente se expõe um 
comentário (alargado) ao artigo em questão. 


Vejam-se portanto os argumentos em que o autor 
se apoia mais decisivamente: 


1.º argumento — «Não se encontra a mínima cor- 
relação estatística entre os valores do coeficiente aero- 
dinâmico C; e os consumos dos carros a 120 Km/h. 
Foi este, com efeito, o resultado duma análise inci- 
dindo sobre uma amostra composta por 15 modelos 
disponíveis no mercado inglês; todos os carros têm 
motores de 1600c.c., os coeficientes C; respectivos 
variam entre 0,38 e 0,47 e os consumos vão desde 7,1 
até 9,0 litros por 100 quilómetros». 


Comentário I — Constatação de facto, com todo o 
aspecto de irrefutável. Torna-se no entanto útil notar 
o critério simplista que parece ter presidido à selecção 
da amostra referida: no elenco dos veículos exami- 
nados acham-se certamente diferenças significativas, 
relativamente a outros factores que não o aerodi- 
nâmico, o que é bem capaz de ter ocultado a corre- 
lação procurada. 

Por outro lado, a simples consideração do consumo 
sem atender ao seu incremento com a velocidade 
(Cio — Coo) mais terá contribuido para tal ocultação. 


2.º argumento — «No campo do aerodinamismo, 
o factor realmente representativo não é o coeficiente 
Ca mas sim o seu produto pela área frontal A do 
veículo. Nesta base, a correlação encontrada mostra 
uma significância relativamente fraca: a amostra em 
questão quando reexaminada sob este novo critério, 
revela que as diferenças de consumo apenas são expli- 
cadas em 45% pelas diferenças no factor Cj. A». 


Comentário II — De novo matéria de facto, mas 
desta vez com mais sólido fundamento que a anterior: 
bem se compreende como a resistência aerodinâmica 
não pode ser indiferente às dimensões do veículo. 
Contudo, sem esquecer o que acima foi sugerido 
quanto aos consumos (incremento mais significativo 
que o valor absoluto) estranha-se que 45% seja con- 
siderado um fraco grau de correlação. Trata-se visi- 
velmente de uma atitude de preconceito que Interessa 
desmontar: na realidade, como os factores interve- 
nientes são 3 ou 4, a conclusão correcta é que o 
aspecto aerodinâmico bem pode ser considerado como 
a componente singular predominante — já que, por 
si só, quase iguala o efeito combinado das restantes. 

Adicionalmente, surge aqui uma verificação digna 
de relevo e que não é sublinhada no texto, a saber: 
quanto maior fôr A maior importância terá o apuro 
da forma para reduzir C;— ou seja, o aspecto forma 
deve merecer mais atenção nos veículos de maior 
porte. 

Extrapolando, pode até afirmar-se que a pouca 
atenção concedida a este aspecto na concepção das 
carroçarias dos autocarros que circulam entre nós 
constitui grave omissão sob uma óptica de conserva- 
ção energética [3]. Possível explicação para o facto 
estará, pelo menos parcialmente, na política de preços 
praticada em Portugal, privilegiando em excesso o 
consumo de gasóleo. 


3.º argumento — «A resistência aerodinâmica, por 
ser proporcional ao quadrado da velocidade, só po- 
derá tomar alguma expressão quando se pratiquem 
consistentemente altas velocidades, o que está longe 
de ser o caso para a maioria dos automobilistas. A 
circulação em auto-estrada não atinge sequer 10% 
da quilometragem total dum condutor britânico mé- 


dio; nos itinerários comuns, também não abundam 
as oportunidades de praticar com segurança níveis 
de velocidade suficientemente elevados.» 


«Nesta base e recordando que o limite legal abso- 
luto está fixado, no Reino Unido, em 112 Kkm/hora 
(70 milhas/h) verifica-se — comparando de novo os 
já referidos 15 veículos, mas agora através dos consu- 
mos a 90 Km horários — que o aspecto aerodinâmico 
apenas explica um terço das diferenças encontradas. 
Na prática, atendendo à diversidade de situações que 
caracterizam o trabalho dum automóvel médio usado 
por um cidadão comum, pode afirmar-se que duma 
redução da resistência aerodinâmica em 1% não re- 
sulta geralmente uma melhoria de consumo superior 
a 0,25%.» 


Comentário HI — Este é de facto o argumento 
decisivo, porventura mais expressivo ainda no caso 
português. A referência aos limites legais, por um 
lado significa que o relevo habitualmente dado ao 
coeficiente aerodinâmico (e também à velocidade 
máxima atingível...) traz implícito um convite ao 
desrespeito pela Lei, que qualquer cidadão cons- 
ciente deve regeitar; por outro, sugere que os nossos 
limites não são afinal dos mais restritos, comparados 
com os que vigoram noutros quadrantes, o que tal- 
vez nos possa dissuadir de esperar maior liberalidade 
neste campo; bem pelo contrário!... 


Tal não será, necessariamente, o caso em países 
menos vulneráveis à factura petrolífera e em que as 
distâncias predominantes imponham maior valoriza- 
ção do factor tempo nas deslocações. 


4.º argumento — «Demonstrada a inutilidade prá- 
tica do coeficiente aerodinâmico, mais significativo 
será preocuparmo-nos com o conhecimento de outros 
factores menos publicitados, nomeadamente o peso, os 
pneus, a transmissão, a cilindrada do motor e o estilo 
de condução. Por exemplo, reduzir o peso em 1% 
pode conduzir a um melhoramento do consumo em 
cerca de 0,7%». 


«Nestas condições, o conselho mais eficaz a trans- 
mitir será o de procurar o carro que tenha o motor 
mais pequeno; menor potência significará menor fluxo 
de combustível e também veículo mais leve.» 


Comentário IV — A conclusão está certa, mas o 
conselho está impreciso. Há muitos modelos ofere- 
cidos no mercado em diferentes versões correspon- 
dendo a motores com diferentes cilindradas sendo 
porém o peso praticamente o mesmo. 

Em tais circunstâncias, o melhor consumo (na 
prática) não corresponde forçosamente à menor cilin- 
drada, embora possa suceder que os consumos con- 
vencionais (segundo as normas da «ECE» — Gene- 
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bra») apontem nesse sentido. A influência dos per- 
cursos ingremes e sinuosos, que entre nós abundam, 
arrastará um apelo mais frequente aos escalões bai- 
xos da caixa de velocidades e aí se perderá rendi- 
mento pelas maiores perdas mecânicas do motor. 
Noutro extremo, as cilindradas demasiado avultadas 
geram igualmente perdas, seja através do convite a 
fortes acelerações e altas velocidades seja pelo fun- 
cionamento a carga parcial nos casos em que se tenha 
de atender a restrições exógenas (p. ex. na circulação 
urbana). 

Interessante é agora notar, a propósito, que, em 
número posterior da revista inglesa acima mencio- 
nada, dá-se notícia da próxima publicação pela 
O.C.D.E. dum «relatório especial» dedicado aos pro- 
blemas do consumo energético dos veículos automó- 
veis e suas perspectivas próximo futuras [4]. 

Ao que se anuncia, ali serão abordadas algumas 
questões candentes como a pressão dos pneus, O 
emprego de materiais plásticos, as transmissões de 
variação contínua, os computadores de bordo e o 


96 


rendimento dos motores face às exigências da con- 
dução automóvel. Sem esquecer os múltiplos factores 
exógenos aos veículos que, na condução urbana, se 
sobrepõem acentuadamente ao efeito das caracterís- 
ticas próprias dos mesmos sobre o seu consumo 
efectivo. | 

Um oportuno documento, sem dúvida, que se 
espera venha a despertar entre os nossos engenheiros 
um maior interesse pelo problema. 
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OPERADORES ASSOCIÁVEIS À RADIAÇÃO 
EM MEIOS EM MOVIMENTO 


M. JOÃO M. MARTINS 


Instituto Superior Técnico (I.S. T.) 
Centro de Electrodinâmica (1. S.T.) 


RESUMO 


Neste artigo ocupamo-nos das equações de onda 
associáveis à radiação em meios sem perdas, homo- 
géneas em movimento uniforme. 

Procura-se reduzir através de transformações ade- 
quadas, c equação à forma usual hiperbólica. Para 
tal definem-se novos operadores e obteve-se uma 
transformação de coordenadas que oferece interesse 
pela sua simplicidade formal. 


ABSTRACT 


In this papel the wave equations that describe 
radiation from sources in lossless moving media are 
studied. We introduce new operators and coordinate 
transformations that reduce the above equations to 
their usual hyperbolic form. 


Considere-se um meio homogéneo, linear sem 
perdas em movimento uniforme em relação a um 
referencial S. As equações que regem os campos são 
as equações de Maxwell. 


oB 
VXE = — —— (1) 
ot 
oD 
VHS (2) 


original recebido para publicação em 20/1/85 


V.D=p (3) 
V.B=0 (4) 
em que p e J são as densidades de carga e corrente. 


Devido ao meio estar em movimento as relações 
constitutivas são as obtidas por Tai (1964) 


D=ec.E+toxH (5) 
B=uva.H-9QXE (6) 


em que « é a diádica 


ds 
o (7) 


o ep 
“ni P (8) 
nº = Le/m€ (9) 
8 (nº—1) 
“ec I-ni (0) 
p=— (11) 
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Supomos que v = ve, e que n$ < 1. 


Estas relações constitutivas são equivalentes às 
obtidas por Minkowski. 


Substituindo-se (5) e (6) em (1) a (4) obtém-se 


VXE=—y da 12 
=-—u 7 (12) 
VaXH = J+e (13) 
Va CB, =P (14) 
VaWH. =0 (15) 
em que 
o q | 
— E Suns jr 16 
Va (7 = ) (16) 
E. =a.E (17) 
H.=e.H (18) 
m=P+02.] (19; 


A equação (15) sugere-nos a possibilidade de es- 
crever 


WH.= "Ya XA (20) 


em que A seria um potencial vector que verifica 
Va HH, = 0. Na realidade prova-se que assim é. 


Por substituição de (20) em (12) obtém-se 


2A 
=—Vab—— 21 


dA 
uma vez que se prova também que VmnX [E + Es = 


JA 
implica — Via P= E + Ea 


Esta relação é análoga à que se obtém para o 
operador usual. 
Recorrendo a (20) e (21) obtém-se de (13) e (14) 


= = Ee ra 
Va x| EV aXA | =-ew TG A- 
= o ob - 
— eW E Va +] (22) 
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tg (0) 2 
m = 0 l/a o (23) 


e existe sempre para n$8 1 


= já 
(VV) + Vad = (24) 


Atendendo à seguinte identidade (Tai) 
= -] 1 = 
v.x( a xA= 2) v) 
(9 .A)= 9.8. VIA | (25) 


que se verifica igualmente para o operador Vm 
e usando 


as equações (22) e (24) tomam a forma 


(Va. %. Va) — EA das ves 
| ot? E 


(26) 
= E | 

(Vm. 4 Vo)AÃ— e pari ro: A)= —u'ai] (27) 

Considerando fontes que se movem segundo z, o 


potencial vector À tem apenas uma componente A, 
e a equação (27) reduz-se a uma equação escalar. 


As equações (26) e (27) revestem então a forma 


- 26 
l nP—pe qd ap 
a a o 4 
ce (I-rb) dt” E 
? JA 
a : (RE SA "a2] 29) 
e d-np o? - Ka. (29) 
em que se substitui 
Pa=+OJ=ap paraJ=Je, (30) 


Podemos pois definir um novo operador 


Pelo que (28) e (29) revestem a forma 


pa 
E im 

E 
[aÃ = —y'8] 


introduzindo a transformação 


x =Vax y=Vãy 27 =2z 
f=t-—or 


podemos exprimir os operadores na forma usual: 


= à 
e(v-0 7) = Wr 
dt 


Õ fi 


E 
eee 
T— 


E 
obtendo-se a equação hiperbólica corrente. 
—u'a?] 


V2 — q? 


nº? =| 


Atendendo à métrica desta transformação 


l 
a £u — Ey = —OC 

E 

Bi A 


Eu = —(1 — co? 
gu=1 
c admitindo a completa covariância das equações de 


Maxwell sob esta transformação é-se conduzido à 
seguinte transformação para os campos 


E' = «(E + co X B) 


d=-1 


B'=aa ..B 
= Q 
D=(D+XH) 
CZ 
H' = aa H 
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VARIAÇÃO DE ENTROPIA 


E REVERSIBILIDADE 


RODRIGO DE ABREU * 


SUMÁRIO 


O «Segundo Principio» da termodinâmica é cons- 
trupdo a partir da verificação de que o aumento de 
energia dum sistema para um dado valor das variá 
veis de deformação aumenta a força que o sistema 
exerce no exterior. Assim sendo o Segundo Princípio 
afirma-se na existência duma relação entre a ener- 
gia, entropia e variáveis de deformação, relação que 
não se postula mas que decorre daquela verificação. 
Desta forma o «Segundo Princípio» constrói-se com 
independência dos chamados Princípios «Zero» e 
«Primeiro». 


SUMMARY 


The «Second Principle» cf Thermodynamics 1s 
built upon the verification that the increase of energy 
in a system for a given value of the deformation 
variables increases the force exerted by the system on 
the outside. This being the case, the Second Prin- 
ciple is affirmed mn the existence of a relation bet 
ween energy, entropy and deformation variables, a 
relation which is not postullated but assumed on the 
basis of the above verification. In this way the Second 
Principle is established independently of the so called 
Zero and First Principles. 


* Instituto Superior Técnico, da UTL e 
Centro de Electrodinâmica, do INIC 
original recebido para publicação em 7/6/85 
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Introdução 


É bem conhecida a dificuldade associada aos 
conceitos de entropia e de reversibilidade. Ainda 
hoje diversos trabalhos têm neste domínio procurado 
evitar contradições, esclarecer conceitos e apurar 
destes quais os essenciais [1, ..., 20]. 

Neste primeiro artigo o conceito de entropia surge 
intimamente ligado à conservação da energia, atra- 
vés dum sistema finito que não mais pode regressar 
ao estado inicial. Isto é a irreversibilidade. No entanto 
a transformação reversível surge com facilidade como 
a transformação limite em que a entropia não varia. 
Nesta fase da exposição diríamos que a reversibili- 
dade ou não reversibilidade está ligada à variação de 
entropia. 

A generalização para situações mais complexas é 
imediata mas neste artigo aflora-se o fundamental. 


O sistema finito que vamos considerar é um gás 
que por simplicidade e para salientar o essencial 
vamos admitir constituído por partículas sem interac- 
ção e cujo número se mantém; admitiremos também 
que nestas condições a pressão do gás aumenta com 
o aumento de energia para um dado volume. 

Admitamos o gás contido num recipiente cujas 
paredes não intervém energeticamente no processo a 
não ser por se deslocar um êmbolo que se encontra 
num dos topos e de massa m.. 

O gás ocupa inicialmente o volume V, à pres- 
soã Pi. 

A pressão Pi é igual à pressão exercida pelo êm- 
bolo e resulta tão somente do facto da massa m, 


estar num campo gravitacional de aceleração g (cons- 
tante). O êmbolo está a uma altura h; (Fig. 1). 

A energia inicial do gás é E, e admitiremos por 
simplicidade e sem perda de generalidade que esta 
é só a cinética: admitiu-se que não havia interacção 
e desprezou-se também a energia associada à exis- 
tência dum campo gravitacional. 

Coloquemos uma massa m;, sobre o êmbolo. 
O sistema evolui e o volume final de equilíbrio é V; 
e a energia do gás é E, Um modo de o fazer é 
suspender m, em contacto e sem pressão por um fio 
que se corta num dado instante. 


FIG. 1 


Admitamos que providencialmente existe uma 
janela J à altura h,. Se assim retirarmos a massa m; 
à altura h, (Fig. 1, 2 e 3) o destacável dá lugar à 
massa m, que se move a altura constante, sem atrito, 
e o êmbolo irá subir para uma altura h, que se de- 
monstrará que não é igual a h;, embora haja a ten- 
tação de admitir que é h, (Ver Apêndice). 

De facto: 


E =Et+meg(h—h)-mg(h/—h) 


Se admitirmos que h; = h; e dado p; = p; ter-se-ia 
que E/= E, o que é absurdo. 

É evidente que dado E/ > E; h, > h;. Isto porque 
se admitiu que a pressão para um mesmo volume 
aumenta com a energia (Fig. 3). 


FIG.2 


Coloquemos à altura h' a massa m” tal que o 
volume final de equilíbrio, após se ter cortado o fio 
de suspensão, seja Vi. 


A energia final E, será 


E =E+mg(y-h)+meg(h'—h) 


Tem-se, portanto, que E; contém a mais que E, 
AE e AÉE' 


AE = m;g (h;—hy) 


sE'=m'g(h!— hj) 


É evidente que sem intervenção dum agente ex- 
terno as massas m, é m' não podem regressar às altu- 
ras iniciais h; e h' respectivamente (Fig. 4) e tão 
pouco o gás pode regressar à energia E,. 


Il — Temos, portanto, que a energia do gás não 
pode ser função apenas do volume. 


ELE(V) 


2 — Verificou-se, também, que não é possível 
fazer regressar o gás ao estado inicial sem 
intervenção duma acção exterior. Portanto a 
transformação é irreversível e quando o gás 
regressar ao volume inicial tem uma energia 
final maior do que a inicial (E, > E;). 


3 — Se a massa m, tender para zero, E,» E. 
Nesta base uma transformação reversível pode 
conceber-se através do acréscimo de massas 
infinitesimais. Para estas transformações 
quando o gás regressar ao volume inicial a 
Energia regressa ao valor iniciál. (É genera- 
lizada a confusão entre transformação rever- 
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sível e quase-estática [1,2,3,6,9, 11, 14,18, 
19, 20]. Tal confusão leva a contradições que 
só desaparecem após a revisão conceptual 
que passa por este esclarecimento [1, 19]. 


Introdução da variável Entropia 


Dado o anteriormente verificado conclui-se que 
o parâmetro V não é um parâmetro suficiente para 
descrever o estado energético. São necessários outros. 
Admitamos que é necessário apenas um. E esse é a 
entropia. Que deste modo e nesta acepção se intro- 
duz aqui. (Admitir a existência dum só parâmetro é 
equivalente à verificação experimental de serem ne- 
cessárias duas variáveis para caracterizar o estado do 
sistema). 


4 — As transformações reversíveis são tais que 


AS = 0 o que é consistente com 3. 


5 — Para as transformações irreversíveis arbitra-se 
AS > 0 para E(S+AS, V) > E(S, V) o que é 
consistente com 2. Feita esta última arbitra- 
riedade é possível afirmar que numa trans- 
formação irreversível a entropia aumenta. 


Generalização para um sistema «infinito» 


Consideremos o elevador representado na Fig. 5. 


FIG.5 


Acrescentada a massa m no andar 1 o elevador 
desce para o andar (i—1). Quando atingir o andar 1, 
através do acréscimo sucessivo de diversas massas m, 
introduz-se na gaveta vazia do andar 1 (Fig. 5) 
uma das massas m e sobe-se para o 2.º andar. Tal é 
possível, sem haver incompatibilidade com a análise 


102 


feita num sistema finito, dado estarmos a conside- 
rar que o sistema tem dimensões tais que as variá- 
veis macroscópicas por unidade de volume não sofrem 
alteração significativa para valores da energia intro- 
duzida no sistema inferior a um dado valor; deste 
modo a descida da massa m de i para 1 não impede 
o regresso do elevador ao 2.º andar, após se ter reti- 
rado m em 1, 

Quando o elevador regressar ao andar 1, através 
de operações sucessivas idênticas à anterior, o resul- 
tado final equivale à descida da massa m de i para 
1 [17], dado todas as gavetas ficarem preenchidas 
excepto a do andar i. O sistema tem, portanto uma 
energia maior do que a inicial. Como não é possível, 
sem intervenção duma acção exterior, fazer com que 
o sistema regresse à energia inicial, nem tão pouco 
fazer com que a massa m regresse ao andar 1 ou 
ainda as duas coisas, o elevador não funciona. (Uma 
forma de aparentemente recuperar as condições ini- 
ciais, mas por isso inútil, é descer as escadas e ir 
buscar a massa m!). 


Conclusão 


Mostrou-se neste trabalho que se considerarmos 
um sistema de energia E e de variável de deforma- 
ção V na presença dum campo gravitacional g, de 
tal modo que a descida duma massa m devido a g 
tenha uma variação de energia potencial gravitacio- 
nal simétrica da que corresponde à subida de m da 
mesma altura, (não existência de irreversibilidade 
associada à posição de m no campo), então a energia 
do sistema não pode ser função apenas de V. De- 
monstrou-se esta afirmação para sistemas que em 
pontos de equilíbrio exerçam uma força crescente 
com a energia para um dado V. Isto é óbvio no 
modelo utilizado recorrendo a conceitos ligados a 
experiências elementares. Esta constatação permite, 
rapidamente, atingir o resultado de que a energia não 
é só função do volume. Se considerarmos transfor- 
mações quaisquer que regressem à configuração ini- 
cial, a energia final é maior do que a inicial [2]. 
Dito doutra forma a transformação é irreversível e 
E(S+AS, V) > E(S, V) sendo S a entropia e podendo 
arbitrar-se AS > 0. 


Ao se admitir a relação anterior entre a força e 
a energia concebe-se a transformação reversível como 
a que se realiza ao longo de pontos de equilíbrio 
muito próximos (quase-estática) de tal forma que o 
sistema possa regressar ao estado inicial. Para tal o 
sistema refere-se à parte do universo para a qual é 
definida a entropia (o exterior é reversível — coa- 
ceito de trabalho associado à subida ou descida de 
pesos) e a transformação reversível é a de variação 
de entropia nula. O 2.º Princípio surge assim da 


relação entre a força e a energia e este é o resul- 
tado fundamental deste trabalho. 


Não se confunde reversibilidade com transforma- 
ção quase-estática nem se introduzem termos como 
adiabático e dissipativo conseguindo-se uma formu- 
lação directa em termos da variável entropia [2, 18]. 

Identificados que sejam os termos calor e ener- 
gia interna [1], para uma mesma configuração a 
mais calor corresponde mais entropia, o que dito 
desta maneira torna evidente a impossibilidade da 
transformação do calor em trabalho regressado que 
seja o sistema à configuração inicial. Tal não impede 
a transformação do calor em trabalho se a variável 
de deformação variar e implica que o trabalho se 
transforma em calor se a variável de deformação 
for a mesma. Na linguagem da entropia (e feita a 


arbitrariedade 
EN 
Vas )y7 


estas últimas afirmações correspondem a que a en- 
tropia dum sistema não pode diminuir. Isto é a 
irreversibilidade. E nesta óptica e adquiridos desta 
forma os conceitos facilmente se compreendem algu- 
mas das dificuldades e contradições exstentes neste 
domínio da física (*) nomeadamente na formulação 
das transformações relativistas das grandezas da ter- 
modinâmica [1], bem assim se adquire uma lingua- 
gem adequada que permite facilmente compreender 
algumas das formulações que ultimamente têm 
aumentado de importância por causa da necessidade 
de optimizar o aproveitamento da energia [20]. Esta 
adequação deve-se ao facto de se conseguir inter- 
pretar duma forma muito directa o segundo Prin- 
cípio da Termodinâmica, evitando-se os conflitos que 
a formulação de Clausius originou e que infelizmente 
tem sido recuperada nas formulações tipo-Caratheo- 
dory [1, 11]. Estas últimas afirmações procuram en- 
caminhar os estudiosos desta matéria no sentido do 
fundamental, eliminando-se as contradições e difi- 
culdades de abordagem conducentes a falsas interpre- 
tações e a questões sem sentido físico [1]. 


(*)«A study of the most general conditions under wich 
heat flow and work have meaning would doubtless be of 
great interest and so for as I am aware has never been 
attempted». — P. W. Bridgman, The Nature of Thermodyna- 
mics (Harper and Broth, New York, 1961). 


Apêndice 


Admitamos que h' < h;. Então existirá uma dada 
massa que se pode retirar a m, e que leva o êmbolo 
a subir até h;. A energia seria E 


E; = E, + m,g(h;—h;) + me g(h;—h9) 
E; = E; + m;g(h;—hy) + me glh;—h)— 
—m/ g(h;—h) 
m' = me — |Am| 
Como m' <m, 
E = E +m;g(h—hy) + A 


A>0 
Logo 
E > E, 


P;< P, para um mesmo V 


O que é absurdo, dado se ter admitido que para 
um mesmo V se E; > E, P; > P.. 
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